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Abstract 

A unified theory of phase transitions and quantum effects in quan- 
tum anharmonic crystals is presented. In its framework, the relationship 
between these two phenomena is analyzed. The theory is based on the 
representation of the model Gibbs states in terms of path measures (Eu- 
clidean Gibbs measures). It covers the case of crystals without translation 
invariance, as well as the case of asymmetric anharmonic potentials. The 
results obtained are compared with those known in the literature. 
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1 Introduction and Setup 

In recent years, there appeared a number of publications describing infucncc of 
quantum effects on phase transitions in quantum anharmonic crystals, where 
the results were obtained by means of path integrals, see [31 El [HI 13 HS1 HH1 
E01 ESI EH E2] ■ Their common point is a statement that the phase transition 
(understood in one or another way) is suppressed if the model parameters obey 
a sufficient condition (more or less explicitely formulated). The existence of 
phase transitions in quantum crystals of certain types was proven earlier, see 
[13 HH H3 HH E0] , also mostly by means of path integral methods. At the same 
time, by now only two works, [6 and |54j . have apperared where both these 
phenomena are studied in one and the same context. In the latter paper, a more 
complete and extended version of the theory of interacting systems of quantum 
anharmonic oscillators based on path integral methods has been elaborated, see 
also [TnillHIH] f° r more recent development, and E3j where the results of El] 
were announced. The aim of the present article is to refine and extend the 
previous results and to develop a unified and more or less complete theory of 
phase transitions and quantum effects in quantum anharmonic crystals, also in 
the light of the results of [531 154] , Note that, in particular, with the help of these 
results we prove here phase transitions in quantum crystals with asymmetric 
anharmonic potentials^, what could hardly be done by other methods. 



This result was announced in |39| . 
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The quantum crystal studied in this article is a system of interacting quan- 
tum anharmonic oscillators indexed by the elements of a crystal lattice L, which 
for simplicity we assume to be a d-dimensional simple cubic lattice H l . The 
quantum anharmonic oscillator is a mathematical model of a quantum particle 
moving in a potential field with possibly multiple minima, which has a suffi- 
cient growth at infinity and hence localizes the particle. Most of the models 
of interacting quantum oscillators are related with solids such as ionic crystals 
containing localized light particles oscillating in the field created by heavy ionic 
complexes, or quantum crystals consisting entirely of such particles. For in- 
stance, a potential field with multiple minima is seen by a helium atom located 
at the center of the crystal cell in bcc helium, see page 11 in [43]. The same 
situation exists in other quantum crystals, He, H2 and to some extent Ne. An 
example of the ionic crystal with localized quantum particles moving in a double- 
well potential field is a KDP-typc ferroelectric with hydrogen bounds, in which 
such particles are protons or deuterons performing one-dimensional oscillations 
along the bounds, see [HI [72l [78l [79]. It is believed that in such substances 
phase transitions are triggered by the ordering of protons. Another relevant 
physical object of this kind is a system of apex oxygen ions in YBaCuO-type 
high-temperature superconductors, see [23 \S7\ [7U [7^ ■ Quantum anharmonic 
oscillators are also used in models describing interaction of vibrating quantum 
particles with a radiation (photon) field, see [321 E21 EH], or strong electron- 
electron correlations caused by the interaction of electrons with vibrating ions, 
see [251 [25] , responsible for such phenomena as superconductivity, charge den- 
sity waves etc. Finally, we mention systems of light atoms, like Li, doped into 
ionic crystals, like KC1. The quantum particles in this system are not necessar- 
ily regularly distributed. For more information on this subject, we refer to the 
survey |34j . 

To be more concrete we assume that our model describes an ionic crystal 
and thus adopt the ferroelectric terminology. In the corresponding physical 
substances, the quantum particles carry electric charge; hence, the displacement 
of the particle from its equilibrium point produces dipole moment. Therefore, 
the main contribution into the two-particle interaction is proportional to the 
product of the displacements of particles and is of long range. According to 
these arguments our model is described by the following formal Hamiltonian 

H=-\Y, J tv-fa<i#)+Y, Ht ' ( L1 ) 

Here the sums run through the lattice L = Z d , d E N, the displacement, qt, of 
the oscillator attached to a given I G L is a ^-dimensional vector. In general, we 
do not assume that the interaction intensities Ju> have finite range. By (•, •) and 
I • I we denote the scalar product and norm in K", M. d . The one-site Hamiltonian 

Ht = Ht + V t (qt) = ^Ptf + faf + Vtfa), a>0, (1.2) 

describes an isolated quantum anharmonic oscillator. Its part H^ r corresponds 
to a (/-dimensional harmonic oscillator of rigidity a. The mass parameter m 
includes Planck's constant, that is, 



m = m ph /h 2 , 



(1.3) 
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where m p h is the physical mass of the particle. Therefore, the commutation 
relation for the components of the momentum and displacement takes the form 

pfqP - qppf = -i5u>&ir > ./•/ 1 "• (1-4) 

For a detailed discussion on how to derive a model like (jl.ip , (|1.2|l from physical 
models of concrete substances, we refer the reader to the survey [72] . 

The theory of phase transitions is one of the most important and spectac- 
ular parts of equilibrium statistical mechanics. For classical lattice models, a 
complete description of the equilibrium thermodynamic properties is given by 
constructing their Gibbs states as probability measures on appropriate config- 
uration spaces. Usually, it is made in the Dobrushin-Lanford-Rucllc (DLR) 
approach which is now well-elaborated, see Georgii's monograph [30j and the 
references therein. In general, the quantum case does not permit such a uni- 
versal description. For some systems with bounded one-site Hamiltonians, e.g., 
quantum spin models, the Gibbs states are defined as positive normalized func- 
tionals on algebras of quasi-local observables obeyng the condition of equilibrium 
between the dynamic and thermodynamic behavior of the model (KMS condi- 
tion), see [19] . However, this algebraic way cannot be applied to the model 
(jl.ip , (|1.2[) since the construction of its dynamics in the whole crystal L is 
beyond the technical possibilities available by this time. In 1975, an approach 
employing path integral methods to describe thermodynamic properties of mod- 
els like (jl.ip , (|1.2j) has been initiated in pQ . Its main idea was to pass from real 
to imaginary values of time, similarly as it was done in Euclidean quantum field 
theory, see [311 167] , and thereby to describe the dynamics of the model in terms 
of stochastic processes. Afterwards, this approach, also called Euclidean, has 
been developed in a number of works. Its latest and most general version is 
presented in [S3} [53] , where the reader can also find an extensive bibliography 
on this subject. The methods developed in these works will be extensively used 
in the present study. 

Phase transitions are very important phenomena in the substances modeled 
by the Hamiltonians (jl.ip . (|1.2p . According to their commonly adopted physi- 
cal interpretation, at low temperatures the oscillations of the particles become 
strongly correlated that produces macroscopic ordering. The mathematical the- 
ory of phase transitions in models like (II. ip , (|1.2|) is based on quantum versions 
of the method of infrared estimates developed in [28]. The first publication 
where the infrared estimates were applied to quantum spin models seems to be 
[21] . After certain modifications this method, combined with path integral tech- 
niques, was applied in [15j [TBI E3 HU [60] to particular versions of our model. 
The main characteristic feature of these versions was a symmetry, broken by 
the phase transition. 

In classical systems, ordering is achieved in competition with thermal fluc- 
tuations only. However, in quantum systems quantum effects play a significant 
disordering role, especially at low temperatures. This role was first discussed 
in 65j. Later on a number of publications dedicated to the study of quan- 
tum effects in such systems had appeared, see e.g., [SS] [82] and the references 
therein. For better understanding, illuminating exactly solvable models of sys- 
tems of interacting quantum anharmonic oscillators were introduced and stud- 
ied, see [B2] [73] Hill EJJ- In these works, the quantity to -1 = fr 2 / m ph was used 
as a parameter describing the rate of quantum effects. Such effects became 
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strong in the small mass limit, which was in agreement with the experimental 
data, e.g., on the isotopic effect in the ferroelectrics with hydrogen bounds, see 
[TS1 175] , see also [57J for the data on the isotopic effect in the YBaCuO-type 
high-temperature superconductors. However, in those works no other quantum 
effects, e.g., those connected with special properties of the anharmonic poten- 
tials, were discussed. At the same time, experimental data, see e.g., the table 
on page 11 in the monograph 18j or the article [75], show that high hydrostatic 
pressure applied to KDP-type ferroelectrics prevents them from ordering. It is 
believed that the pressure shortens the hydrogen bounds and thereby changes 
the anharmonic potential. This makes the tunneling motion of the quantum 
particles more intensive, which is equivalent to diminishing the particle mass. 
In [5J [HI |5] > a theory of such quantum effects in the model (|1.1[) . (|1.2[) . which ex- 
plains both mentioned mechanisms, was buit up. Its main conclusion is that the 
quantum dynamical properties, which depend on the mass to, the interaction 
intensities Jw, and the anharmonic potentials Vi, can be such that the model 
is stable with respect to phase transitions at all temperatures. 

As was mentioned above, the aim of this article is to present a unified descrip- 
tion of phase transitions and quantum stabilization in the model (|1.1|) . Q1.2p . 
mostly by means of methods developed in [53] [SJ . We also give here complete 
proofs of a number of statements announced in our previous publications. The 
article is organized as follows. In Section[2j we briefly describe those elements of 
the theory developed in [531 HI] which we then apply in the subsequent sections. 
In Section[3] we present the theory of phase transitions in the model (jl.l[) , (|1.2[) . 
We begin by introducing three definitions of a phase transition in this model 
and study the relationships between them. Then we develop a version of the 
method of infrared estimates adapted to our model, which is more transpatrent 
and appropriate than the one employed in [6]. Afterwards, we obtain a suffi- 
cient conditions for the phase transitions to occur in a number of versions of 
the model (jl.ip . Q1.2p . This includes also the case of asymmetric anharmonic 
potentials Vt which was never studied before. At the end of the section we 
make some comments on the results obtained and compare them with similar 
results known in the literature. Section [4] is dedicated to the study of quantum 
stabilization, which we understand as the suppression of phase transitions by 
quantum effects. Here we discuss the problem of stability of quantum crystals 
and the ways of its description. In particular, we introduce a parameter (quan- 
tum rigidity), responsible for the stability and prove a number of statements 
about its properties. Then we show that under the stability condition which we 
introduce here the correlations decay 'in a proper way', that means the absence 
of phase transitions. The relationship between the quantum stabilization and 
phase transitions are also analyzed. In the simplest case, where the model is 
translation invariant, scalar (y — 1), and with the interaction of nearest neigh- 
bor type, this relation looks as follows. The key parameter is 8dm Jd^, where 
d is the lattice dimension, J > is the interaction intensity, and > is 
determined by the anharmonic potential V (the steeper is V the smaller is 
Then the quantum stabilization condition (respectively, the phase transition 
condition) is 8dm J^l < 1, see (]4.32jl . (respectively, 8dm J$l > <f>(d), see (13.71)) 
and (|4.33p ). Here 4> is a function, such that <fi(d) > 1 and cf>(d) — ► 1 as d — > +oo. 
We conclude the section by commenting the results obtained therein. 
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2 Euclidean Gibbs States 

The main element of the Euclidean approach is the description of the equilibrium 
thermodynamic properties of the model , (II. 2p by means of Euclidean Gibbs 
states, which are probability measures on certain configuration spaces. In this 
section, we briefly describe the main elements of this approach which are then 
used in the subsequent parts of the article. For more details, we refer to |54j . 

2.1 Local Gibbs states 

Let us begin by specifying the properties of the model described by the Hamil- 
tonian The general assumptions regarding the interaction intensities Jw 

are 

Ju> = Jt'i > 0, J«=0, J = supV^ Jw < oo. (2.1) 

i J, 

In view of the first of these properties the model is ferroelectric. Regarding 
the anharmonic potentials we assume that each Ve : M 1 ' — > R is a continuous 
function, which obeys 

A v \x\ 2r +B v < Vi(x) < V(x), (2.2) 

with a continuous function V and constants r > 1, Ay > 0, By G K- In certain 
cases, we shall include an external field term in the form 

V e (x) = V e °(x) - (h, x), h e R v , (2.3) 

where V® is an appropriate function. 

Definition 2.1 The model is translation invariant if Vi = V for all £, and the 

interaction intensities Jw are invariant under the translations o/L. The model 
is rotation invariant if for every orthogonal transformation U G 0(v) and every 
I, Vi(Ux) — Vi(x). The interaction has finite range if there exists R > such 
that Jw = whenever \£ — i'\ > R. 

If Vi = for all £, one gets a quantum harmonic crystal. It is stable if Jo < a, 
see Remark l2.15l below. 

By A we denote subsets of the lattice L; we write A (e L if A is non-void 
and finite. For such A, by |A| we denote its cardinality. A sequence of subsets 
A <s L is called cofinal if it is ordered by inclusion and exhausts the lattice L. If 
we say that something holds for all £, we mean it holds for all £ € L; sums like 

mean X^gl- W e also use the notations K + = [0, +oo) and No = NU {0}, N 
being the set of positive integers. 

Given AgL, the local Hamiltonian of the model is 

Ha = -\ J "' ' fa* ««') + E H * ( 2 - 4 ) 

which by the assumptions made above is a self-adjoint and lower bounded op- 
erator in the physical Hilbert space L 2 (R ly l A l). For every (3 = 1/ksT, T being 
absolute temperature, the local Gibbs state in A <g L is 



qa(A) = tr&ce[Aexp(-PH A )]/Z A , A e £ A , 



(2.5) 
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where 

Z A = trace[exp(-/3ff A )] < oo (2.6) 

is the partition function, and Ca is the algebra of all bounded linear operators 
on L 2 (R I/ ' A I). Note that adjective local will always stand for a property related 
with a certain AgL, whereas global will characterize the whole infinite system. 

The dynamics of the subsystem located in A is described by the time auto- 
morphisms 

= exp(itH\)Aexp(—itH\), (2-7) 

where t G R is time. Given n G N and A\, . . . ,A n G (£a> the corresponding 
Green function is 

G%,... An (ti, ...,*») = (?A K (Ar) • • • afj4„)] , (2.8) 

which is a complex valued function on R n . Each such a function can be looked 
upon, see as the restriction of a function G\ A analytic in the domain 

V% = {( Zl , . . . , z n ) G C" | < Q( Zl ) <■■■< %{z n ) < (3}, (2.9) 

and continuous on its closure. The corresponding statement is known as the 
multiple-time analyticity theorem, see [HE], as well as |41] for a more general 
consideration. For every n G N, the subset 

{(*!, . . . , z n ) G V n I ft( Zl ) =■■■ = di(z n ) = 0} (2.10) 

is an inner uniqueness set for functions analytic in M, see pages 101 and 352 
in [66] . This means that two such functions which coincide on this set should 
coincide everywhere on DS. 

For a bounded continuous function F : R"I A I — > C, the corresponding multi- 
plication operator F G £a acts as follows 

(FV)(sb) =F(x)i/>(x), v e L 2 (R v|a| ). 

Let 5a C £a be the set of all such operators. One can prove (the density 
theorem, see [SUES]) that the linear span of the products 

at 1 (F 1 )---at(F n ), 

with all possible choices of n G N, ti, . . . , t n G R, and . . . , F n G # A i is dense 
in Ca in the cr-weak topology in which the state (|2.5p is continuous as a linear 
functional. Thus, the latter is determined by the set of Green functions Gp p 
with n G N and Fx, ... ,F n G 3a- The restriction of the Green functions ^ 
to the imaginary-time sets (|2.10j) are called Matsubara functions. For 

Tl < t 2 < • • • < r „ < /3, (2.11) 

they are 

r^,... ;F „(ri, . . . ,r n ) = G^...^^, . . . ,tr n ). (2.12) 

Since (I2.10|l is an inner uniqueness set, the collection of the Matsubara functions 
(|2.12|) with all possible choices of n G N and Fi , . . . , F n G 3a determines the 
state (|2.5|) . The extensions of the functions (|2.12p to [0, f3] n are defined as 

^F lt ...,F n ( T li • ••! r n) = r F<r (n) ( T CT(l) ) ■ ■ • J T er(n)) 
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where a is the permutation such that T^-m < T^-p) " • • T o-(n)- One can show 
that for every 9 £ [0,0], 

r F 1 ,...,F„(n + A, • • • ,r n + 6) = T^. Fn ( Tl , r n ), (2.13) 
where addtion is modulo 

2.2 Path spaces 

By (j2.8[) . the Matsubara function (|2.12p can be written as 

^F 1 ,...^(n,...,r n )= (2.14) 



trace 



P e -(T 2 -Tl)H A p ' g-(T3-T 2 )i?A . ^ g-(r„ + l-T n )H A 



/^A, 



where r„ + i — (3 + t\ and the arguments obey (12. lip . This expression can be 
rewritten in an integral form 

^F 1 ,..,F n (r 1 ,...,T n )= I F 1 {tu A {T 1 ))---F n (u JA {T n ))v A {dLu A ), (2.15) 

that is the main point of the Euclidean approach. Here is A is a probability 
measure on the path space Q A which we introduce now. The main single-site 
path space is the space of continuous periodic paths (temperature loops) 

C l3 = {4>e C([0, -> W) | 0(0) = <t>(0)}. (2-16) 

It is a Banach space with the usual sup-norm || • \\cp- F° r an appropriate <f> G Cp, 
we set 

KM)=F- sup „><,, (2.17) 

T,T'e[0,/3] t^t' l r — r 1/3 

where 

|t - r'|/3 = min{|r - t'|;/3 - |t - r'|} (2.18) 

is the periodic distance on the circle Sp ~ [0,/3]. Then the set of Holder- 
continuous periodic functions, 

C£ = {</> G C/3 | K a {4>) < (2.19) 

can be equipped with the norm 

H\\o S = \<l>(0)\+K tr (<f,), (2.20) 

which turns it into a Banach space. Along with the spaces Cp, CZ, we shall 
use the Hilbert space L\ = L 2 {Sp — > K",dT), equipped with the inner product 
(•, -)i2 and norm || • \\ L 2 . By B(Cp), B(L^) we denote the corresponding Borel a- 

algebras. In a standard way, see page 21 of 59J and the corresponding discussion 
in [S3], it follows that 

Cp G B{Ll) and B(Cp) = B(L^) n Qj. (2.21) 

Given A C L, we set 

Ha = {^a = {ijJt)teA I ^ G Cp}, (2.22) 
= f2 L = = (^)^ e jL | ^ € C^}. 
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These path spaces are equipped with the product topology and with the Borcl 
cr-algebras B(Q.\). Thereby, each Cl\ is a complete separable metric space, 
called Polish space, its elements are called configurations in A. For A C A', the 
juxtaposition lja> — uja X u>a'\A defines an embedding fl\ <—> Qa' by identifying 
uja € f^A with uja x Oa'\a G f^A'- By P(Oa), V(d) we denote the sets of all 
probability measures on (Ua, B(fl\)), (tt,B(£l)) respectively. 



2.3 Local Euclidean Gibbs measures 

Now we construct the measure va which appears in (12.15[) . A single harmonic 
oscillator is described by the Hamiltonian, c.f., (|1.2j) . 

fl »"-s:g(^) > + 5 w "- <2 ' 23) 

It is a self-adjoint operator in the space L 2 (R y ), the properties of which are well- 
known. The operator semigroup exp(— riJ^ 1 ar ), t E Sp, defines a /3-periodic 
Markov process, see |42j . In quantum statistical mechanics, it first appeared 
in R. H0egh-Krohn's paper [35 . The canonical realization of this process on 
(Cp,B(C/3)) is described by the path measure which can be introduced as fol- 
lows. In the space L|, we define the following self-adjoint Laplace-Beltrami type 
operator 

d 2 



A = I -m-^ + a 1 <g> I, (2.24) 

where I is the identity operator in M 1 ". Its spectrum consists of the eigenvalues 

A; = m(2irl/l3) 2 + a, I e Z. (2.25) 

Therefore, the inverse A -1 is a trace-class operator on and the Fourier 
transform 



exp 

r 2 



%($,<!>) IA xW) = exp{ -\{A- l i),tl}) Ll \ (2.26) 



defines a zero mean Gaussian measure \ on (^^(L^)). Employing the eigen- 
values (|2.25p one can show that, for any p G N, 



,')f'x^) < 0^±f> (| 



\u{T)-u,{S)\* X {to)< (£:) 'Ir-r^. (2.27) 



Therefrom, by Kolmogorov's lemma (see page 43 of [HHj) it follows that 

X(C|) = 1, for all a e (0,1/2). (2.28) 

Thereby, x(C^) = 1; hence, with the help of (|2.2ip we redefine x as a measure 
on (Cp,B(Cf3)), possessing the property (|2.28p . We shall call it H0egh-Krohn's 
measure. An account of the properties of x can be found in [BJ. Here we 
present the following two of them. The first property is obtained directly from 
Fernique's theorem (see Theorem 1.3.24 in [22]). 
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Proposition 2.2 (Fernique) For every a £ (0,1/2), there exists X a > 0, 
which can be estimated explicitely, such that 



exp 



(M0IIo ? ) X(<ty) < oo. (2.29) 



The second property follows from the estimate (|2.27|) by the Garsia-Rodemich- 
Rumsey lemma, see [33]. For fixed a £ (0, 1/2), we set 

s*M= sup | K | T) ~y 2 L *e (0,0/2), 

r,r': 0<|t— r'| a <i9 [ l T — T l/3 J 

(2.30) 

One can show that, for each a and it can be extended to a measurable map 
E,j : C/3 — * [0, +oo]. 

Proposition 2.3 (Garsia-Rodemich-Rumsey estimate) Given a £ (0,1/2), 
let p £ N be such that (p — l)/2p > a. Then 

S» x (d W ) < D(a,p,v)m- p ^ l - 2 °\ (2.31) 

where m is the mass M.3\) and 

2W"(1 + 1/^ 2?T>/2 + l) 
D{CT ' P ' U) = (p-l-2ap)(p-2ap) ' I>/2) ' ^ 

The H0egh-Krohn measure is the local Euclidean Gibbs measure for a single 
harmonic oscillator. The measure va £ V(£Ia), which is the Euclidean Gibbs 
measure corresponding to the system of interacting anharmonic oscillators lo- 
cated in A (e L, is defined by means of the Feynman-Kac formula as a Gibbs 
modification 

Va((^a) = exp [-Ja(wa)] XA(dwA)/^A (2.33) 

of the 'free measure' 

xA(d^ A )-n^ d ^)- ( 2 - 34 ) 

^SA 

Here 



1 V 



(2.35) 



is the energy functional which describes the interaction of the paths u>j>, £ £ A. 
The normalizing factor 



N A = / exp[-/ A (wA)]XA(dw A ) (2.36) 

is the relative partition function, whereas the Feynman-Kac representation of 
the partition function (|2.6p is 

Z A = N A Z\ aI , (2.37) 
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where 



ghar dcf trace exp 



eeA 



exp 



-03/2)Vo/m 



„|A| 



1 — exp I —(3^Ja/m 



Now let us summarize the connections between the description of the sub- 
system located in A (g L in terms of the states (|2.5[) and of the Euclidean Gibbs 
measures (|2.33[) . By the density theorem, the state qa is fully determined by the 
Green functions (|2.8[) corresponding to all choices of n £ N and F\, . . . , F n £ 3a- 
Then the multiple-time analyticity theorem leads us from the Green functions 
to the Matsubara functions (|2.12j) , which then arc represented as integrals over 
path spaces with respect to the local Euclidean Gibbs measures, see (|2.15|) . 
On the other hand, these integrals taken for all possible choices of bounded 
continuous functions F%, . . . , F n fully determine the measure v\. Thereby, we 
have a one-to-one correspondence between the local Gibbs states (|2.5j) and the 
states on the algebras of bounded continuous functions determined by the local 
Euclidean Gibbs measures (|2.33p . Our next aim is to extend this approach to 
the global states. To this end we make more precise the definition of the path 
spaces in infinite A, e.g., in A = L. 



2.4 Tempered configurations 

To describe the global thermodynamic properties we need the conditional dis- 
tributions 7TA(do;|£), A <s L. For models with infinite-range interactions, the 
construction of such distributions is a nontrivial problem, which can be solved 
by imposing a priori restrictions on the configurations defining the corresponding 
conditions. In this and in the subsequent subsections, we present the construc- 
tion of such distributions performed [54 . 

The distributions 7TA(doj|£) are defined by means of the energy functionals 
I\(uj\^) describing the interaction of the configuration u> with the configuration 
£, fixed outside of A. Given A <s= L, such a functional is 

JaMO = Ja(wa) - Yl Ju\ojUi')l%. u£tl, (2.38) 

eeA, i'eA-= 

where I a is given by (|2.35|) . Recall that u> — uia x loa"] hence, 

JaMO = 7 A (w a x A .|0 A x f A „). (2.39) 

The second term in (I2.38|) makes sense for all £ £ fi only if the interaction 
has finite range, see Definition 12.11 Otherwise, one has to impose appropriate 
restrictions on the configurations £, such that, for all I and uj £ f2, 

Ju> ■ \(o>t,St')i%\ < oo. (2.40) 

t' 

These restrictions are formulated by means of special mappings (weights), which 
define the scale of growth of {||& ||l2 }^£l- Their choice depends on the asymp- 
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totic properties of Jw, \l — i'\ — > +00, see (I2.1[) . If for a certain a > 0, 

sup V JweKp(a\£-£' |) < 00, (2.41) 

then the weights {w Q (^, <")} a£ i are chosen as 

w a (£,£') =exp(-a\£-e'\), I = (0,a), (2.42) 

where a is the supremum of a > 0, for which (|2.41|1 holds. If the latter condition 
does not hold for any a > 0, we assume that 



sup 



Y^'hv ■{l + \£-t'\) ad , (2.43) 



for a certain a > 1 . Then we set a to be the supremum of a > 1 obeying (|2.43j) 
and 

w a (e i £') = (l+e\£-£'\)- ad , (2.44) 

where e > is a technical parameter. In the sequel, we restrict ourselves to 
these two kinds of J«/. For more details on this item, we refer the reader to 

Given a€l and u> £ f2, we set 



1/2 

^|M| 2 l2Wq (0,£) 



(2.45) 



and 

n a = {lo e n I |M| Q < 00}. (2.46) 

Thereby, we endow fl a with the metric 

^ 1 + 11^-^11^ 

which turns it into a Polish space. The set of tempered configurations is defined 
to be 

0* = P| fl a . (2.48) 

We endow it with the projective limit topology, which turns it into a Polish 
space as well. For every a 6 X, the embeddings f2* <— > 51 Q f2 are continuous; 
hence, f2„, fi* £ 2?(r2) and the Borel a-algebras B(fl a ), Bffi) coincide with the 
ones induced on them by B(fl). 

2.5 Local Gibbs specification 

Let us turn to the functional ()2.38|) . By standard methods, one proves that, for 
every a € X, the map fl a x tt a > I\(to\^) is continuous. Furthermore, for any 
ball B a (R) = {cj e fi Q I p Q (0,w) < i?}, i? > 0, one has 

inf 7aM£) > —00, sup |7a(u;|£)| < +00. 
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Therefore, for AgL and £ € fl*, the conditional relative partition function 

N A (£)= [ exp[-/ A K x A c|£)] XA (dc; A ) (2.49) 

is continuous in £. Furthermore, for any i? > and q£I, 

inf 2V A (£) > 0. 
ieB a (R) 

For such £ and A, and for B £ B(tt), we set 

*a(B\0 = — i-r / exp[-/ A ( WA x 0a=|£)]IbK x £ A c) XA (du;A), (2.50) 

where Is stands for the indicator of B. We also set 

7T A (-|0 = 0, for £ € n\n*. (2.51) 
From these definitions one readily derives a consistency property 

7r A (B|w)7rA'(dw|0 = 7r A '(S|0» AcA', (2.52) 

which holds for all B € 6(0) and 

The local Gibbs specification is the family {7r A } A gL. Each 7r A is a measure 
kernel, which means that, for a fixed £ G O, 7r(-|£) is a measure on (fi, B(f2)), 
which is a probability measure whenever £ £ f2*. For any S £ $(Q), 7r A (S|-) is 
£?(il)-measurable. 

By Cb(O a ) (respectively, Cb(fi*)) we denote the Banach spaces of all bounded 
continuous functions / : £l a — > K (respectively, / : O* — * R) equipped with the 
supremum norm. For every a £ X, one has a natural embedding Cb(f2 Q ) °- ► 
C b (O*) . Given a € X, by W Q we denote the usual weak topology on the set of 
all probability measures V{^l a ) defined by means of Cb(f2 a ). By VV* we denote 
the weak topology on ^(O*). With these topologies the sets V(fl a ) and Vffi 1 ) 
become Polish spaces (Theorem 6.5, page 46 of [55]). 

By standard methods one proves the following, see Lemma 2.10 in [S"4"] , 

Proposition 2.4 (Feller Property) For every a e X, A <g L, and any / £ 

Cb(f2 Q ), the function 

n a 3^7r A (f\0 (2.53) 

def 1 Z 1 

= 7777V / /( WA x ^ Ac ) ex P [- j a( w a x M£)] XA(dw A ), 
A'aU) Jn A 

belongs to Cb(f2 Q ). 27ie linear operator f i— ► tta(/|-) *s a contraction on Cb(Q a )- 
Note that by (|23U|) . for (eSl'.ael, and / € C b (O a ), 

*a(/|0 = / /(w)7r A (da;|0. (2.54) 



Recall that the particular cases of our model were specified by Definition 12.11 
For B E B(n) and U £ 0{v), we set 

Uu = (ULu e ) eeh UB = {Uuj | uj e B}. 
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Furthermore, for a given £q, we set 

t lo (uj) = (w £ _ £o ) £eL , U {B) = {tt (w) \ uj £ B}. 
Then if the model possesses the corresponding symmetry, one has 

ir A (UB\US) = 7r A (B|0, n h+t {t t {B)\t t ^)) = tt a (B|£), ( 2 -55) 
which ought to hold for all U, £, B, and £. 

2.6 Tempered Euclidean Gibbs measures 

Definition 2.5 .A measure ji € is called a tempered Euclidean Gibbs mea- 

sure if it satisfies the Dobrushin-Lanford-Ruelle (equilibrium) equation 

[ Tr A (B\uj)fi{doj) = n(B), for all A(eL and B e 0(0). (2.56) 
Jo 

By C7* we denote the set of all tempered Euclidean Gibbs measures of our model 
existing at a given f3. The elements of Q t are supported by 0*. Indeed, by (|2.50p 
and (T2~5T|) tt a (0 \ O'|0 = for every A <e L and £ € fi. Then by p3S)l . 

j u(Q\f2 t ) = 0. (2.57) 

Furthermore, 

At ({w S fi* | W S L : w £ S C£}) = 1, (2.58) 

which follows from (|2.28|) . (|2.29|) . If the model is translation and/or rotation in- 
variant, then, for every U € 0{y) and £ 6 L, the corresponding transformations 
preserve 5*. That is, for any fj, £ Q t 1 

Qu{n) d ^ f ^oU^eG", 6 t {n) not^eQ*. (2.59) 

In particular, if 5* is a singleton, its unique element should be invariant in the 
same sense as the model. From Proposition 12.41 one readily gets the following 
important fact. 

Proposition 2.6 For each a £ X, every W Q - accumulation point fi £ 'P(O t ) of 
the family {7r A (-|£) | A gL, £ S f2*} is an element of Q t . 

Now let us pay some attention to the case where the model (|1.1[) . p.2|) is 
translation invariant. Recall that the lattice L = "L d is considered as an additive 
group. For £o £ L, A g L, and cj G O, we set 

A +4 = + U € A}; t <0 (w) = (£|°)*el, C' = (2.60) 

Furthermore, for B G 0(0), we set 

t t {B) = {t t (u) | w G B}. (2.61) 

Clearly, ^(B) G B(fl) and t^(f2*) = 0* for all & 

Definition 2.7 A probability measure /i £ V(ft) is said to be translation in- 
variant if for every I and B £ B(£l), one has ii(ti(B)) = j-i{B). 
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As was mentioned above, the Gibbs specification {ttaIahl of the translation 
invariant model is translation invariant, that is, it has the property (|2.55[) . 

Remark 2.8 The translation invariance of the Gibbs specification does not 
mean that each probability kernel 7ta as a measure is translation invariant. 
Moreover, it does not mean that all the Euclidean Gibbs measures defined by 
this specification are translation invariant. One can only claim that if the set 
£/* consists of one element only, this element ought to translation invariant. 

Set 

B mv = { Be g(Q) | W • tt ( B } = B } ; (2.62) 

which is the set of all translation invariant events. By construction, f2* 6 B lnv . 
We say that fieV(Cl) is trivial on B inv if for every B € B inv , one has fi(B) = 
or n{B) = 1. By 7 3inv (f2) we denote the set of translation invariant probability 
measures on (Q,B). 

Definition 2.9 A probability measure fi £ "P mv (0) is said to be ergodic (with 
respect to the group h) if it is trivial on B lnv (Q). 

Ergodic measures are characterized by a mixing property, which we formulate 
here according to see Theorem III. 1.8 on page 244. For L £ N, we set 

Al = (-L, L] d n Z rf , (2.63) 

which is called a box. For a measure [i and an appropriate function /, we write 

0% = J /dM (2-64) 

Proposition 2.10 (Von Neumann Ergodic Theorem) Given /i S V mv (Cl), 
the following statements are equivalent: 
(i) /i is ergodic; 

(li) for all f,g e L 2 (tt,ii), 

Ao \KT\ { ^ {L fWoitfaMM - (/>/i ' (.9)m) |=0. (2.65) 

Proposition 2.11 If the model is translation invariant and Q is a singleton, 
its unique element is ergodic. 

Now we give a number of statements describing the properties of Q t . More 
details can be found in [5"4"| . 

Proposition 2.12 For every (3 > 0, the set of tempered Euclidean Gibbs mea- 
sures Q t is non-void, convex, and W'- compact. 

Recall that the Holder norm || • \\c° was defined by (|2.20p . 

Proposition 2.13 For every a € (0,1/2) and x > Q, there exists a positive 
constant C such that, for any t and for all fi € Q l , 

[ exp(\ a \\Lu e \\ 2 c „ +x\\u t \\l 3 ) h(6lj) < C, (2.66) 

where \ a is the same as in h2.29\) . 



2 EUCLIDEAN GIBBS STATES 



16 



In view of (|2.66[) . the one-site projections of each fj, £ Q l are sub-Gaussian. The 
constant C does not depend on t and is the same for all [i £ Q t , though it may 
depend on a and h. The estimate (|2.66p plays a crucial role in the theory of 
the set Q t . 

According to [30] certain Gibbs states correspond to the thermodynamic 
phases of the underlying physical system. Thus, in our context multiple phases 
exist only if Q t has more than one element for appropriate values of (3 and the 
model parameters. On the other hand, a priori one cannot exclude that this set 
always has multiple elements, which would make it useless for describing phase 
transitions. The next statement which we present her^l clarifies the situation. 
Let us decompose 

V t = V u + V 2/ , (2.67) 
where Vi,t S C 2 {W) is such that 

- a < b = f inf inf (^(x)y, y) l\y\ 2 < oo. (2.68) 

As for the second term, we set 

< <S d = sup ( sup V 2 t(x) - inf V 2 e(x)\ < oo. (2.69) 

Its role is to produce multiple minima of the potential energy responsible for 
eventual phase transitions. Clearly, the decomposition (|2.67|) is not unique; its 
optimal realizations for certain types of Vi are discussed in section 6 of [13] . 
Recall that the interaction parameter Jo was defined in (|2.ip . 

Proposition 2.14 The set Q t is a singleton if 

e 136 < (a + b)/j . (2.70) 

Remark 2.15 The latter condition surely holds at all (3 if 

(5 = and J < a + b. (2.71) 

If the oscillators are harmonic, (5 = 6 = 0, which yields the stability condition 

J < a. (2.72) 

The condition \2. 70\ ) does not contain the particle mass m; hence, the property 
stated holds also in the quasi-classical limits m — * +oo. 

By the end of this subsection we consider the scalar case v = 1. Let us introduce 
the following order on Q l . As the components of the configurations lo S fl are 
continuous functions uji : Sp — > ffi", one can set lo < lo if loi(t) < Cogir) for all I 
and t. Thereby, 

K+ffl*) = {/ e C b (n*) | f(oj) < f(u>), if lo < Co}, (2.73) 
which is a cone of bounded continuous functions. 

2 C.f., Theorem 3.4 in [54], Theorem 2.1 in [H], and Theorem 4.1 in [13] . 
3 More details on this limit can be found in \E\. 
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Proposition 2.16 If for given n,fi& Q t , one has 

(/>m =</>£= fOT a11 f&K + (Q% (2.74) 

then fj, = fl. 

This fact allows for introducing the FKG-order. 
Definition 2.17 For /J.,fj,G Q l , we say that [i<fi, if 

(/)m<</)m= fOT a11 /€«+(«*)• (2.75) 

Proposition 2.18 TTie setQ 1 possesses maximal /i+ and minimal [i^ elements 
in the sense of Definition \2.17\ These elements are extreme; they also are 
translation invariant if the model is translation invariant. If Vi(—x) = Vt(x) 
for all I, then ft+{B) = fi-(-B) for all B G B(fi). 

The proof of this statement follows from the fact that, for / 6 K + (il t ) and any 
AdL, 

< </WlC). whenever £ < (2.76) 
which one obtaines by the FKG inequality, see [51) . By means of this inequality, 
one also proves the following 

Proposition 2.19 The family {7TA(-|0)}AgL has only one W -accumulation point, 
/j,q, which is an element of ' G t ■ 

2.7 Periodic Euclidean Gibbs measures 

If the model is translation invariant, there should exist <j> '■ Nq — > R + such that 
Ju> =<f>(\ti \e d -e d \). (2.77) 

For the box |(2j55)| . we set 

JA d J: f ^x-^U,...,!^-^), (2.78) 

where 

|4 - = minfl^ - 4| ; L - & - ^|}, j = 1, . . . ,d. (2.79) 
For £, ^' £ A, we introduce the periodic distance 

\e-e'\A = ^\h-£' 1 \l + --- + \£ d -e d \ 2 L . (2.80) 

With this distance the box A turns into a torus, which one can obtained by 
imposing periodic conditions on its boundaries. Now we set, c.f., (|2.35p . 




(2.81) 



and thereby, c.f., (|2.33D . 

^r(d^A) = ex P [-/r^A)]xA(dwA)/^r r J (2.82) 

NF* = I cxp[-/r>A)]XA(da, A ). 
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By means of (|2.78p we introduce the periodic Hamiltonian 

and the corresponding periodic local Gibbs state 

g p A CI (A) = trace[Aexp(-/3 J ffP cr )]/trace[exp(-/3i7P G1 ')], A e £ A . (2.84) 

The relationship between the measure v^ cr and this state is the same as in the 
case of i/ A and g\. 
Set, c.f., (12301) . 

= i / ex P [-T r K)] Ib(wa x OaOxa^a), (2.85) 
which is a probability measure on O*. Then 

7rP or (d( WA x WA „)) = ^"(d^A) J] ^, (da*), (2.86) 

feA c 

where 0^' is the zero element of the Banach space Cp. Note that the projection 
of 7r^ er onto 17a is ^a"- 

Let £box be the sequence of all boxes (|2.63p . Arguments similar to those 
used in the proof of Lemma 4.4 in [51] yield the following 

Lemma 2.20 For every a £ 2 and a S (0, 1/2), there exists a constant C > 
such that, for all boxes A, 

J nt fe IMIc^(0,^ < cr (d^) < C. (2.87) 

Thereby, the family {^^}KeC bm is W* -relatively compact. 

Let be the family of W'-accumulation points of {T^'lAeCbox- 

Proposition 2.21 It follows that A4 C £7* . The elements of M., called periodic 
Euclidean Gibbs measures, are translation invariant. 

The proof of this statement is similar to the proof of Proposition ^. 61 It can be 
done by demonstrating that each \i S M. solves the DLR equation (|2.56[) . To 
this end, for chosen A gL, one picks the box A containing this A, and shows 
that 

/^A(-|£K° r (d£)^M-), if in W'. 

Here both convergence are taken along a subsequence of £box- 

2.8 The pressure 

In the translation invariant case, one can introduce a thermodynamic function, 
which contains important information about the thermodynamic properties of 
the model. This is the pressure, which in our case up to a factor coincides with 
the free energy density. As our special attention will be given to the dependence 
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of the pressure on the external field h, c.f. (|2.3[) . we indicate this dependence 
explicitely. For A d L, we set, see (|2.49[) . 

P\(h, £) = ±- log N A (h,£), eefi*. (2.88) 

To simplify notations we write p A {h) = P\{h, 0). Thereby, for ^ £ 5', we set 

p£(h)= [ PA(h,Z)n{d$). (2.89) 



Furthermore, we set 

p A eT (h) = ±-logNr i (h). (2.90) 
If, for a cofinal sequence £, the limit 

p»(h) ^ limp^h), (2.91) 
exists, we call it pressure in the state /x. We shall also consider 

p(h) d £ limp A (h), p^(h) H f lim pl c \h). (2.92) 

C £box 

Given I = (h, . . .l d ), I' = . . ,l' d ) e L = Z d , such that lj < I'j for all j = 
1, . . . , d, we set 

r = {£ G L | lj < £j < l' v for all j = 1, . . . , d}. (2.93) 

For this parallelepiped, let &(T) be the family of all pair- wise disjoint translates 
of T which cover L. Then for A <g L, we let 7V_(A|r) (respectively, iV+(A|r)) 
be the number of the elements of <S(T) which are contained in A (respectively, 
which have non-void intersections with A). 

Definition 2.22 A cofinal sequence C is a van Hove sequence if for every T, 
(a) limAT_(A|r) = +oo; (6) lim (iV_(A|r)/iV + (A|r)) = 1. (2.94) 

One observes that £box is a van Hove sequence. It is known, see Theorem 3.10 
in [54], that 

Proposition 2.23 For every h € K and any van Hove sequence C, it follows 
that the limits h2.91}) and \2. 92]) exist, do not depend on the particular choice 
of C, and are equal, that is p(h) — p pcl '(h) — p^(h) for each fi € Q l . 

Let the model be rotation invariant, see Definition 12.11 Then the pressure 
depends on the norm of the vector h £ HI". Therefore, without loss of generality 
one can choose the external field to be (h, 0, . . . , 0), h £ K. For the measure 
(|2.33|) . by vi ' we denote its version with h = 0. Then 



/ ex p ( h Yl [ 



■Y.\(/>) - .V A (ll) / exp j /,• > ] lu [ P{t)At\ V f\dw K ). (2.95) 



2 EUCLIDEAN GIBBS STATES 



20 



The same representation can also be written for iV? er (/i). One can show that 
the pressures p A (h) and p A or (/i), as functions of h, are analytic in a subset of C, 
which contains M. Thus, one can compute the derivatives and obtain 

^p A (h) = (3M A (h), ^PT(h) = PMrW, (2-96) 

where 

M A (h) d = f -^jErf], Mr\h) d ^ (2.97) 

are local polarizations, corresponding to the zero and periodic boundary condi- 
tions respectively. Furthermore, 

J^PaW ( 2 -98) 



2 

1 



2|A| 



W ( W «(r)-£>«(r))dr 



^A(dWA)^A(dw A ) > 0. 



The same can be said about the second derivative of p A CI (/i). Therefore, both 
p A (h) and p^ cr (/i) are convex functions. For the reader convenience, we present 
here the corresponding properties of convex functions following 69] , pages 34 - 
37. 

For a function p : R — > K, by y±(i) we denote its one-side derivatives at a 
given t £ K. By af most countable set we mean the set which is void, finite, or 
countable. 

Proposition 2.24 For a convex function <p : R — > R, it follows that: 

(a) the derivatives <p'±(t) exist for every t € R; 

the set {t £ R | tp' + (t) ^ p>'_(t)} is at most countable; 

(b) for every t£R and 9 > 0, 

¥>'_(*) < ^ ¥>'_(* + 0) < + (2-99) 

(c^ £/ie point-wise limit ip of a sequence of convex functions {</? ra }neN 
is a convex function; if p and all ip n 's are differentiable at a 
given t, (p' n (t) — > <p'(t) as n — > +oo. 

Proposition 2.25 The pressure p{h) , see Provosition \ 2.2$[ is a convex func- 
tion o//tel. Therefore, the set 

Tl d = {h E K I p'_(h) < p'+(h)} (2.100) 

is at most countable. For any h €E 1Z C and any van Hove sequence C, it follows 
that 

limM A (h) = lim M^ cr (/i) = pr x p'{h) = f M{h). (2.101) 

£ £box 

By this statement, for any h € 1Z C , the limiting periodic state is unique. In the 
scalar case, one can tell more on this item. The following result is a consequence 
of Propositions |2~2"51 and [2~TH1 
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Proposition 2.26 If v = 1 and p(h) is differentiable at a given h GM., then G t 
is a singleton at this h. 

Returning to the general case v G N we note that by Proposition ^. 251 the global 
polarization M(h) is a nondecreasing function of h G TZ C ; it is continuous on 
each open connected component of 1Z C . That is, M(h) is continuous on the 
intervals (a_,a+) C 1Z C , where a± are two consecutive elements of 1Z. At each 
such o±, the global magnetization is discontinuous. One observes however that 
the set TZ C may have empty interior; hence, M(h) may be nowhere continuous. 

In the sequel, to study phase transitions in the model with the anharmonic 
potentials V of general type, we use the regularity of the temperature loops 
and Proposition 12.31 Let the model be just translation invariant, i.e., the 
anharmonic potential has the form (|2.3p , where V° is independent of £. Let us 
consider the following measure on Cp: 



A(dw) = 



1 



■ exp 



N, 



exp 



V Q (uj(T))dT X (du), 



(2.102) 



where x is H0egh-Krohn's measure. For a box A, we introduce the following 
functions on Q A 



Y A (u; A ) = 
Then from (|2.90[) one gets 

P r\h) = \ og N 
i 



2 E J «' E / 

E 



j = 1, . . . ,ZA 



eeA 



r)dr, 



(2.103) 




exp 



II A(dw^k04) 



eeA 



As the measure (|2.102p is a perturbation of the H0egh-Krohn measure, we can 
study the regularity of the associated stochastic process by means of Proposition 
531 Fix some p£N\ {1} and a G (0, 1/2 - l/2p). Thereby, for <d G (0, /?), one 
obtains 



C/3 



EP(cu)X(dLo) < e 



By being as in 



where, see (|2.32l) . 



By Proposition 1 2 . 31 this yields 
{%)x<D v {a, u,p)m- p ^ X - 2a \ 



(2.105) 



, def 2W 1 )(l + l/ap) 2 f 



(p — 1 — 2pcr) (p — 2per) 



2Pexp(-/3B v )r(^/2+p) 
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For c > and n G N, n > 2, we set 

C^njc) = {w e C/j ±w (j) (fc/?/n) > c, j = 1,. ..,v; k = 0,1,... n}. (2.106) 

For every n G N, ji,...,j n G {1, ■ • • > an d Ti,.. . ,t„ € [0, /?], the joint distri- 
bution of w( Ji )(ti), . . . , w^' n '(T n ) induced by H0egh-Krohn's measure x 1S Gaus- 
sian. Therefore, x((7 (n;c)) > 0. Clearly, the same property has the measure 
(|2.102jl . Thus, we have 

E(n;c) d - min{A(C + (n;c)) ;A(C - (n;c))} > 0. (2.107) 
For e 6 (0, c), we set 

A(c; e) = G C/3 | S^w) < (c - £) 2 (/3/n)" 2CT }, (2.108) 

B ± (e,c) = A(c;e)p|C* ± (n;c). 

Then for any r € [0, /?], one finds k G N such that |r — kj3/n\ < (3/n, and hence, 
for any j = 1, . . . , u, 

\^(r)~^(kf3/n)\ < [E 0/n (uj)] 1/2 (f3/nr, 

which yields ±lo^\t) > e if lo G B ± {e,c). Let us estimate A[B ± (e, c)]. By 
(|2.105jl and Chebyshev's inequality, one gets 



Air'A-iu-:.-.) < n2 J {c _ £)2p (^ /n )x 

— [mn(c — e) 2 ] p 



Thereby, 



A[B ± (e,c)] = A[C ± (n;c)\(C /3 \A( C ;e))] (2.109) 



E(n;c)-A(C^\A(c; £ )) 



> 

> E(n;c) 



[mn(c — e) 
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def , . 
= 7N, 



which is positive, see (|2. 107)) . for all 



dof P { D v (a,v,p) \ 1/P 

m - m * = ^F^A s(»;<0 J ■ ( } 

This result will be used for estimating the integrals in (|2.104p . 



3 Phase Transitions 

There exist several approaches to describe phase transitions. Their common 
point is that the macroscopic equilibrium properties of a statistical mechanical 
model can be different at the same values of the model parameters. That is, 
one speaks about the possibility for the multiple states to exist rather than the 
transition (as a process) between these states or between their uniqueness and 
multiplicity. 
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3.1 Phase transitions and order parameters 

We begin by introducing the main notion of this section. 

Definition 3.1 The model described by the Hamiltonians il.l]) , il.ty) has a 

phase transition if \Q | > 1 at certain values of (3 and the model parameters. 

Note that here we demand the existence of multiple tempered Euclidean Gibbs 
measures. For models with finite range interactions, there may exist Euclidean 
Gibbs measures, which are not tempered. Such measures should not be taken 
into account. Another observation is that in Definition 13.11 we do not assume 
any symmetry of the model, the translation invariance including. If the model 
is rotation invariant (symmetric for v — 1, sec Definition 12. ip . the unique el- 
ement of Q t should have the same symmetry. If l^/'l > 1, the symmetry can 
be 'distributed' among the elements of G t - In this case, the phase transition is 
connected with a symmetry breaking. In the sequel, we consider mostly phase 
transitions of this type. However, in subsection 13 . 51 we study the case where the 
anharmonic potentials Ve have no symmetry and hence there is no symmetry 
breaking connected with the phase transition. 

If the model is translation invariant, the multiplicity of its Euclidean Gibbs 
states is equivalent to the existence of non-ergodic elements of G t , see Corollary 
12.111 Thus, to prove that the model has a phase transition it is enough to 
show that there exists an element of which fails to obey (|2.65[) . In the case 
where the model is not translation invariant, we employ a comparison method, 
based on correlation inequalities. Its main idea is that the model has a phase 
transition if the translation invariant model with which we compare it has a 
phase transition. 

Let us consider first the translation and rotation invariant case. Given I and 
j = 1, . . . , i/, we set 



The right-hand side in (|3.1|) does not depend on r due to the property (|2.13|) . 
To introduce the Fourier transformation in the box A we employ the conjugate 
set A* (Brillouin zone), consisting of the vectors p — (p%, . . . ,Pd), such that 




(3.1) 




, . . . , d. 



(3.2) 



Then the Fourier transformation is 




(3.3) 




In order that u>f (r) be real, the Fourier coefficients should satisfy 



4^(r)=^(r). 
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By the rotation invariance of the state (•) !y p<= r , as well as by its invariance with 
respect to the translations of the torus A, it follows that 



Thus, we set 



U DDI 



E ^ 



Xp,e'-i) 



(3.4) 



(3.5) 



'eA 

= JL V D A 
IAI ^ p 



pGA, 



One observes that D A can be extended to all p G (— 7r,7r] d . Furthermore, 



and 



1 1 peA» 



£'GA 

' ' pGA, 



For M A = (u^ eA , <E 



(ma,-D A «a) 12(a) = E D u' u t u t- 



(3.6) 



(3.7) 



(3.8) 



O'GA 



- E 



E«. f^\r)d, 



> 0. 



Thereby, the operator Z) : Z 2 (A) — > ^ 2 (A) is strictly positive; hence, all its 
eigenvalues D A are also strictly positive. 

Suppose now that we are given a continuous function B : (— it, ir] d — > (0, +oo] 
with the following properties: 



(i) / B{p)Ap < oo, 

J(-Tr,7r] d 

(ii) D£<B(p), for all A, \ {0}, 
holding for all boxes A. Then we set 

1 



(3.9) 



But 



(2tt) c 



(-7r,7r] d 



B(p)cos(p,£-l')dp, £,£' eL, (3.10) 



and 



(3.11) 



peA,\{o} 



We also set Bf t = if cither of £, f belongs to A c . 

Proposition 3.2 For every 1,1' , it follows that Bffi — > as i — ► +oo. 
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Proof: By (|3.9[) . B(p) cos(p,£ — £') is an absolutely integrable function in the 
sense of improper Riemann integral. The right-hand side of (|3.11[) is its integral 
sum; thereby, the convergence stated is obtained in a standard way. □ 

From claim (i) of (13. 9j) by the Riemann-Lebesgue lemma, see page 116 in 
[55j . one obtains 

lim B w = 0. (3.12) 

|£-£'|-»+oo 

Lemma 3.3 For every box A and any £,(.' S A, it follows that 

D&,>(Db-B£)+B&,. (3.13) 
Proof: By (J32J), (|3~TTj) . and claim (ii) of one has 

Dh-D&, = E ^sin 2 (p,£-/) 

^ TXT E %)sin 2 (p,£-0 



l A l 

H A R A 



which yields (l3~T3| . □ 

For e 5', we set, c.f., (l3~Tj) . 



D£,=/3/ ((^(r),^(r'))) M dr'. (3.14) 



o 



Corollary 3.4 For every "periodic fi £ Q t , it follows that 

D&,>(Dfr-Bu)+B u >, (3.15) 

holding for any I, £' . 

Proof: For periodic fx £ Q t , one finds the sequence {L,i} ri eN C N, such that 
7r^ r =>■ /x as ri — * +oo, see Proposition ^. 211 This fact alone does not yet mean 

that Df£ n — > Dfc,, what we would like to get. To prove the latter convergence 
one employs Lemma 12.201 and proceeds as in the proof of claim (b) of Lemma 
5.2 in [54 . Then (pHS]) follows from ([3~T3) and Proposition OH □ 

One observes that the first summand in (|3.15jl is independent of I, whereas 
the second one may be neither positive nor summable. Suppose now that there 
exists a positive i? such that, for any box A, 

r4 > 0. (3.16) 

Then, in view of (|3.12|) . the phase transition occurs if 

•d > B u . (3.17) 

For certain versions of our model, we find the function B obeying the conditions 
(|3"T5|1 and the bound (|3~TT|) . Note that under (|3~T75)) and (|3TT7)l by HHg) it follows 
that 

A TO E ^ - Jft„ ^ E ^ > 0- (3-18) 
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Let us consider now another possibilities to define phase transitions in transla- 
tion invariant versions of our model. For a box A, see l|2.63p . we introduce 



Pa = 



1 



'eA 



(/3|A|)2 



and set 



P = limsupPA £ 

L — >+oo 



(3.19) 



(3.20) 



Definition 3.5 The above P is called the order parameter. If P > for given 
values of (3 and the model parameters, then there exists a long range order. 

By standard arguments one proves the following 

Proposition 3.6 // iTO^j) and ([X77p hold, then P > 0. 

The appearance of the long range order, which in a more 'physical' context is 
identified with a phase transition, does not imply the phase transition in the 
sense of Definition 13.11 At the same time, Definition 13.11 describes models with- 
out translation invariance. On the other hand, Definition 13.51 is based upon 
the local states only and hence can be formulated without employing Q t . Yet 
another 'physical' approach to phase transitions in translation invariant models 
like Ijl.lj) . ()1.2|) is based on the properties of the pressure p(h), which by Propo- 
sition 12.231 exists and is the same in every state. It does not employ the set Q t 
and is based on the continuity of the global polarization (|2.101[) . that is, on the 
differentiability of p(h). 

Definition 3.7 (Landau Classification) The model has a first order phase 
transition if p'(h) is discontinuous at a certain h*. The model has a second 
order phase transition if there exists h* G W such that p'(h) is continuous but 
p"(h) is discontinuous at h — h*. 

Remark 3.8 Like in Definition ] 3. 11 here we do not assume any symmetry of 
the model (except for the translation invariance). As p(h) is convex, p'(h) is 
increasing; hence, p"(h) > 0. The discontinuity of the latter mentioned in 
Definition \3. 7| includes the case p"(h*) = +oo, where the polarization M(h) at 
h = /i* grows infinitely fast, but still is continuous. 

The relationship between the first order phase transition and the long range 
order is established with the help of the following result, the proof of which can 
be done by a slight modification of the arguments used in |24j . see Theorem 
1.1 and its corollaries. Let {/i n }AfeN (respectively, {M n } n ^) be a sequence 
of probability measures on K (respectively, positive real numbers, such that 
limM„ = +oo). We also suppose that, for any jel, 

f(y)= Urn -Llog />V„(d M ) (3.21) 



exists and is finite. As the function / is convex, it has one-sided derivatives 
/±(0), see Proposition 
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Proposition 3.9 (Griffiths) Let the sequence of measures {fJ- n }NGfi be as above. 
If /+(0) = /-(0) — 4 1 (is., f is differentiable at y = §), then 



lim 

n — > + oo 



g(u/M n )n n (du) = g(4>), 



(3.22) 



for any continuous g : R — * R, such that \g(u)\ < Ae^'"' with certain A, x > 0. 
Furthermore, for each such a function g, 



n — >+oo 



limsup / g(u/M n )fj, n (du) < 



ze[/i(o),/' (o)] 



<?(*)• 



In particular, if f'_(0) = —fL{0), then for any k GN, 

/|(0) > limsup ( [(u/M n ) 2k Hn(du)) 

Write, c.£, (|2~95)l . 

Nl cl \h) = Nl cr 



l/2fc 



r (0)/ expffcW' 

^ V teA " 



c^(r)dr]^ POr (d^), 



.O.pcr, 



(3.23) 



(3.24) 



(3.25) 



where i/°' per is the local periodic Euclidean Gibbs measure with h — 0. Now let 
{L„}„ g n c N be the sequence such that the sequences of local measures {i/^ cl } 
and {z^T } converge to the corresponding periodic Euclidean Gibbs measures 
li° and /i respectively. Set 



3u e 



(r)dr = u 



(3.26) 



Clearly, each such X n is measurable and isomorphic to R. Let n € N, be 

7 A 



the projection of {^°' pCI } on to this X n . Then 



p{h)= P (0) + f(h), 



(3.27) 



where / is given by (f3~2T1) with such ,u„ and M„ = |A L J = {2L n ) d . Thereby, 
we apply (|3.24p with k — 2 and obtain 



P+(0) > /3 lim sup JP Ai „- 

n — >+oc * 

Thus, in the case where the model is just rotation and translation invariant, the 
existence of the long range order implies the first order phase transition. 

Consider now the second order phase transitions in the rotation invariant 
case. For a € [0, 1], we set, c.f., (|3.19D . 



P 



(a) 



|A| 



l+a 



£ / ^( r ) dr 



v^(duj A ), 



(3.28) 



where A is a box. Then P[ = Pa and, as we just have shown, the existence of 
a positive limit (|3.20p yields a first order phase transition. 
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Proposition 3.10 If there exists a G (0,1), such that for a sequence {L n }, 
there exists a finite limit 

lim P[ a) d = > 0. (3.29) 

Then the model has at h = a second order phase transition. 
Proof: We observe that 



Pi a) = vp'L(0)/P 2 \A\ a . 
Then there exists c > 0, such that 

Pl Ln {0) > c|A L „| a , for all n G N. 
As each p"^ is continuous, one finds the sequence {<5 n } rl gN such that 5 n J, and 

Pl Ln {h) > -Ak Ln \ a , for all h G [0, 5 n ] and neN. (3.30) 
If p"{0) were finite, see Remark 13.81 one would get 

p"(Q)= lim U («5 n )-^ (0)1 /<5„, 

n— *+oo L T1 n J 

which contradicts (|3.30jl . □ 

Proposition 13.101 remains true if one replaces in (|3. 28[1 the periodic local 
measure ^ er by the one corresponding to the zero boundary condition, i.e., 
by i>a- Then the limit in (|3.29[) can be taken along any van Hove sequence 
C. We remind that Proposition 13.101 describes the rotation invariant case. The 
existence of a positive P^ with a > may be interpreted as follows. According 
to the central limit theorem for independent identically distributed random 
variables, for our model with J«< = and V& = V, the only possibility to 
have a finite positive limit in (|3.29[) is to set a = 0. If P^ < oo for nonzero 
interaction, one can say that the dependence between the temperature loops is 
weak; it holds for small Jo. Of course, in this case P^ = for any a > 0. If 
p(°0 g e ts positive for a certain a G (0,1), one says that a strong dependence 
between the loops appears. In this case, the central limit theorem holds with 
an abnormal normalization. However, this dependence is not so strong to make 
p' discontinuous, which occurs for a = 1, where a new law of large numbers 
comes to power. In statistical physics, the point at which P^ > for a G 
(0, 1) is called a critical point. The quantity P^ is called susceptibility, it gets 
discontinuous at the critical point. Its singularity at this point is connected with 
the value of a for which P^ > 0. The above analysis opens the possibility to 
extend the notion of the critical point to the models which are not translation 
invariant. 

Definition 3.11 The rotation invariant model has a critical point if there exist 
a van Hove sequence C and a G (0, 1) such that 



lim-—— — / / loAt)At 



v A (duj A ) > (3.31) 



at certain values of the model parameters, including h and (3. 

Note that by Proposition ^. 101 it follow that in the translation invariant case the 
notions of the critical point and of the second order phase transition coincide. 
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3.2 Infrared bound 

Here, for the translation and rotation invariant version of our model, we find 
the function B obeying ()3.9p . 

For a box A, let E be the set of all unordered pairs (£, £'), i, £' £ A, such that 
\l — £'\a — 1) see (|2.80l) . Suppose also that the interaction intensities (|2.78[) are 
such that J w = J > if and only if (£, £') £ E and hence the measure (I2.82p 
can be written 

t^° r (dw A ) = — i— exp J -- IK ~ ^'llz,| o-A(dw A ), (3.32) 

M ' \ {t,t')EE J 

where 

cr A (dtJA) (3.33) 
= exp UdJ2 IKIIil ~J2 [ ^M T )) dr ) XA(dw A ), 

and 

Y A (0)= / exp (-- ^ I <r A (dw A ). (3.34) 

With every edge (£, £') £ E we associate bw £ L 2 and consider 

Y\(b) = / exp I ^ - w# - b u >\\ 2 L 2 J cr a (dc^A ) • (3.35) 

By standard arguments, see [46j and the references therein, one proves the 
following 

Lemma 3.12 (Gaussian Domination) For every b = (bw)(i,i')eE, bw £ 
Lg, it follows that 

Y A (b)<Y A (0). (3.36) 
Let Xe be the real Hilbert space 

X E = {b= (bw)(e,e>)eE I bw £ L 2 g }, (3.37) 

with scalar product 

{b,c) XE = ^ {bw,cw)Ly (3-38) 
(i,£')eE 

To simplify notations we write e = (£, £'). A bounded linear operator Q : Xe — > 
Xe may be defined by means of its kernel Q 3 e 3 e , (t, t'), j,j' = l,...,v, e, e' G £7, 
and t,t' £ [0,(3]. That is 



(3.39) 
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Let us study the operator with the following kernel 



(3.40) 



where the expectation is taken with respect to the measure p.32[) . This operator 
in positive. Indeed, 



(b, Qb)x E = 



^2 uj e >,b u >) L 2 



> 0. 



The kernel ()3.40j) can be expressed in terms of the Matsubara functions; thus, 
as a function of r, t' , it has the property (|2.13p . We employ the latter by 
introducing yet another Fourier transformation. Set 



K = {k = (2tt//3)k I k £ Z}, 



e fc(r) = 

The transformation we need is 



(3- 1 / 2 coskr, iffc>0; 
-p-Wsinkr, if < 0; 
y/2/0, if k = 0. 



4 j) (r)e fe (r)dr, 



(3.41) 
(3.42) 

(3.43) 
(3.44) 



Then the property (j2~T5|) yields, c.f., (|3~4"|) 

(cDf (fc)<$V)>tf~ =0 if Mfc' or j ^ f. 



Taking this into account we employ in (|3.40p the transformation (|3.43[) and 
obtain 



keic 



with 



(hA)^A)( k ) = { |4f W-^f (*)J ' \^t'{k)-^>[k) 



■So) I 



■Si) i 



(3.46) 



In view of the periodic conditions imposed on the boundaries of the box A the 
latter kernel, as well as the one given by (|3.40p , are invariant with respect to 
the translations of the corresponding torus. This allows us to 'diagonalize' the 
kernel (13.46|) by means of a spatial Fourier transformation (|3.2|) . (|3.3p . Then 
the spacial periodicity of the state (■} v ,v" yields 



(w^(p,A)(a^(p' > jfe)) I/ p- =0 if p + pV0. 



(3.47) 
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Taking this into account we obtain 



}{i lA ){i 2A )(k) = X^fo^C-P.fc))^ (3-48) 
peA, 

ApA) _ e »(p,*i)^ /|A| 1/2 
' e -»Wa) _ e «(-p,4) N ] /|A| 1/2 . 



Since the summand corresponding to p — equals zero, the sum can be re- 
stricted to A* \ {0}. This representation however cannot serve as a spectral 
decomposition similar to (|3.45p because the eigenfunctions here are not normal- 
ized. Indeed, 

(e<^ - e<^'A /IAI 1 / 2 x (e-^ - e"^)) = 2£(p) 

where 

d 

£(p) d = f £[1-co8 Pj ]. (3.49) 

3=1 

Then we set 

o-«,(p) = (e 1 ^ - e^ £ '>) /VPW, P eA,\ {0}, (3.50) 

and 

Q(p,fc) = 2£(p)<<2^(p,A0fl W) (-p,fc)>^> P eA,\ {0}. (3.51) 
Thereby, 

Q(e 1 r 1 ){e 2 ,£' 2 )(T,T / ) = (3.52) 
= X! X! Q ( yP^) a iit' 1 < yP) (J i2i' 2 {~P) e k{T)e k {T'), 

p£A,\{0} fcGK 

which is the spectral decomposition of the operator (|3.39|) . Now we show that 
the eigenvalues (|3.5ip have a specific upper bouncffl 

Lemma 3.13 For every p£ A, \ {0} and k G K,, the eigenvalues VS.51]) obey 
the estimate 

Q(p,k)<l/J, (3.53) 
where J is the same as in \3. 32\) . From this estimate one gets 

{^X Pl k)^X- P ,k))vr ^ 2JIM' peAA {0} - (3 ' 54) 

Proof: The estimate in question will be obtained from the Gaussian domination 
(|3.36p . For ( G 1 and a given b £ Xe, we consider the function <f>{t) = Y\(tb). 
By Lemma T3. 121 <fi"(0) < 0. Computing the derivative from (|3.35[) we get 



<f/ , (Q) = J(b,Qb) Xa -\\b\\% 



E ' 



4 Their natural lower bound is zero as the operator (13.391 1 is positive 
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where the operator Q is defined by its kernel (|3.40[) . Then the estimate (|3.53p 
is immediate. □ 

By ([33]), (JS3H), and ([3~5T]) . we readily obtain 

((& p (t),u)- p (t'))) v p<* = Z. ( , E Qfe 1 ' fc ) cos [ fc ( T ~ T % P7^0, 



which yields, see (13.51) and (|3.53p , 

G $ = m^°^Jkry p ^°- (3 ' 55) 

Comparing this estimate with (|3.9p we have the following 

Corollary 3.14 If the model is translation and rotation invariant with the near- 
est neighbor interaction, then the infrared estimate \3.9\) holds with 

5(P)= 2^)' PeH^U }' B(0)=+oo. (3.56) 

3.3 Phase transition in the translation and rotation in- 
variant model 

In this subsection, we consider the model described by Corollarv l3.14l First we 
obtain the lower bounds for 

{(Ui{T),U)l{T))) v v*r, 

from which we then obtain the bounds (|3.16p . In the case where the anharmonic 
potential has the form 

V(u) = -b\u\ 2 + b 2 \u\ 4 , b>a/2, 6 2 >0, (3.57) 

a being the same as in (|1.1|> . the bound (|3. 16|) can be found explicitly. We begin 
by considering this special case. 

Lemma 3.15 Let V be as in ^3.57\ ). Then, for every A g L, 

((Mr),Mr))) vr >^f^ =»*■ (3-58) 

Proof: Let A be a self-adjoint operator, such that the expressions below make 
sense. Then 

Ql er ([A,[Hr\A}}) (3.59) 
= q p ™ a {AHl cr A + AHl cr A - AAHl cr - H? CI AA) 

= E i4»'i 2 (^r-^r r ){exp[-/3i?n-exp[-/3£;r]} 

> 0. 

Here E per , s £ N are the eigenvalues of the periodic Hamiltonian (|2.83p . A ss i 
are the corresponding matrix elements of A, and g^ r is the periodic local Gibbs 
state (|2.84p . By the Euclidean representation, 



(Mr),^)))^ = J2 ( (^(o)) ) = E^a 



w\ 2 ' 
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Then we take in (|3.59p A = pp , j — 1, . . . ,v, make use of the commutation 
relation (|1.4[) . take into account the rotation invariance, and arrive at 



fi?([A,[HF,A]]) 



_ per 



eTX H 6 + a + 2b ^ + 4fo 2(<zf) 2 J (3-60) 

2 ' 

-2b + a + Ab 2 (v + 2) 



> 0, 

which yields 1(335]) , □ 

Now we consider the case where V is more general as to compare with (|3.57[) . 



Lemma 3.16 Let the model be translation and rotation invariant, with nearest 
neighbor interaction. Then, for every 9 > 0, there exist positive to, and J*, 
which may depend on (3, 9, and on the potential V , such that, for to > to* and 
J > J*, 

((^(r),^(T)))„ r >6. (3.61) 
Proof: Let us rewrite (|2.104j) 

pT{J) = \o g N p 

+ pql°g(^ exp[y A ( WA )]nA(d^)j, (3.62) 



where we indicate the dependence of the pressure on the interaction intensity 
and have set h = since the potential V should be rotation invariant. Clearly, 
p A cl (J) is convex; its derivative can be computed from (|3.62[) . Then we get 

rn ((pwhj,)^ = JjjpT( j ) ( 3 - 63 ) 

> pT(J)-pT(0) 



where E is the same as in (I3.37|) . By the translation invariance and (|2.13[) . one 
gets 

((we,ue>) L *) u p«r < ({{ut,ut)Li) v ^ + <(^',^')l=)„p») / 2 
= (N,^)^)^ =^((^(r),^(r))) !/ per. 

Then we choose e, c, and n as in (|2. 109|) . apply this estimate in (|3.63p . and 
obtain 

0Jd{{ut{T),wt{T)))^ > ^logj [ exp[y A (a; A )]TTA(d^)l 

> /3Jisde 2 + vlog~f(m). (3.64) 
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For m > to* given by (|2.110p , 7(to) > and the latter estimate makes sense. 
Given 8 > 0, one picks e > y/0/v and then finds J* such that the right-hand 
side of the latter estimate equals 9 for J = J* . □ 

To convert (|3.58p and (|3.61|) into the bound (|3.16p we need the function 
/ : [0, +oo) — * [0, 1) defined implicitly by 

/(tttanhu) = u^ 1 tanhu, for u > 0; and /(0) = 1. (3.65) 

It is differentiable, convex, monotone decreasing on (0, +oo), such that tf(t) — > 
1. For t > 6, /(t) s» 1/i to five-place accuracy, see Theorem A. 2 in [23]. By 
direct calculation, 

/'(mt) 1 t-u(1-t 2 ) , 

-77 — r = ; — r, o\, r = tanhtt. 3.66) 

J(UT) UT T + U(l — T z ) 



Proposition 3.17 For every fixed a > 0, the function 

<f>(t) = taf{t/a), t > (3.67) 
is differentiable and monotone increasing to a 2 as t — > +oo . 
Proof: By ([336]) . 

, 2ar{\ - t 2 ) , 

q> (t) = — > 0, ut — wtanhu = t/a. 

t + u(l — T Z j 

The limit a 2 is obtained from the corresponding asymptotic property of /. □ 

Next, we need the following fact, known as Inequality of Bruch and Falk, see 
Theorem IV.7.5 on page 392 of [H] or Theorem 3.1 in [33]. 

Proposition 3.18 Let A be as in H3.59]) . Let also 

b{A) = p- 1 / g p A cr {Aexp[-T J ffP or ]Aexp[r J ffP or ]} dr, 
g(A) = g r (A 2 ) ; c(A) = tf? {[A, [[3H P a\ A}}} , 



Then 



btA) (3.68) 



where f is the same as in \3.65\) . 
Set 

J(d) - T^vJ / ^T-y (3-69) 

where £(jp) is given by (|3.49|) . The exact value of J(S) can be expressed in terms 
of complete elliptic integrals, see [53] and also [37] for more recent developments. 
For our aims, it is enough to have the following property, see Theorem 5.1 in 
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Proposition 3.19 For d > 4, one has 

d^ <J ^<j A c7(d) < d L -r ^ 

where a(d) — > 1/2 as d — > +oo. 

Recall that to is the reduced particle mass (| 1 . 3|) . 

Theorem 3.20 Lei <i > 3, the interaction be of nearest neighbor type, and the 
anharmonic potential be of the form \3.51\ ), which defines the parameter $* . Let 
also the following condition be satisfied 

Sm'&lJ > J(d). (3.71) 
Then for every (3 > /3» , where the latter is the unique solution of the equation 

2£Jtf*/(/?/4mtf*) = J{d), (3.72) 
the model has a phase transition in the sense of Definition 
Proof: One observes that 

[q?\lHr,q^}} = l/m, I e A. (3.73) 
Then we take in (|3.68|) A = qr^ and obtain 



b(A) >((4 3) (0))\ r /- 



4m<(^' 3 (0)) >^ r 



By Proposition ^. 171 d f ([3 / Amd) is an increasing function of Thus, by (|3.58p 
and |3~T|) . 

£>« > f3 2 vd*f{(3/Am$*), (3.74) 

which yields the bound (|3.16p . Thereby, the condition (i) in (|3.17p takes the 
form 

t?*/ (/3/4mi?*) > J{d)/2/3J. (3.75) 
By Proposition 13. 17[ the function 

<t>{[3) = 2/3Ji?*/(/3/4mtf*) 

is monotone increasing and hits the level J(d) at certain /3*. For (3 > /3*, the 
estimate (|3.75p holds, which yields |C/g| > 1. □ 

One observes that /(/3/4TO-)?*) — ► 1 as to — * +oo. In this limit, the condi- 
tion (|3.7ip turns into the corresponding condition for a classical model of <f) 4 
anharmonic oscillators, Now let us turn to the general case. 

Theorem 3.21 Let d > 3, the interaction be of nearest neighbor type, and the 
anharmonic potential be rotation invariant. Then, for every f3 > 0, there exist 
to* and J* > 0, which may depend on [3 and on the anharmonic potential, such 
that \Q t \ > 1 for m > to* and J > J*. 

Proof: Given positive [3 and 9, the estimate (|3.6ip for big enough m and J. 
Then one applies Proposition l3.181 which yields that the condition (i) in (|3.17|) 
is satisfied if 

9f{f3/Am6) > J(d)/2f3J. 

Then one sets to* to be as in (|2.110p and J* to be the smallest value of J for 
which both (|3.6ip and the latter inequality hold. □ 
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3.4 Phase transition in the symmetric scalar models 

In the case v = 1, we can extend the above results to the models without 
translation invariance and with much more general Jw and Vg. However, certain 
assumptions beyond (|2.ip and (|2.2p should be made. Suppose also that the 
interaction between the nearest neighbors is uniformly nonzero, i.e., 

inf J U ' d = J>0. (3.76) 

\i-i'\=i 

Next we suppose that all Vgs are even continuous functions and the upper bound 
in (|2.79p can be chosen to obey the following conditions: 

(a) for every £, 

V(ug) — Vg{ug) < V(ug) — Ve(ug), whenever uj < u'j; (3.77) 

(b) the function V has the form 

r 

V( Ui ) = b {s) uf s ; 2bM < -a; 6 (s) > 0, s > 2, (3.78) 

s=l 

where a is as in fll.lj) and r > 2 is either positive integer or infinite; 

(c) if r = +oo, the series 

^) = E 2 J( S _ lv 6(s)r " 1 ' ( 3 - 79 ) 

s=2 ^ '■ 

converges at some d > 0. 
Since 26*- 1 -* + a < 0, the equation 

a + 26 (1) + <?(??)= 0, (3.80) 

has a unique solution > 0. By the above assumptions, all Vg are 'uniformly 
double-welled'. IfVg(ug) — vg(uj) and vg are differentiable, the condition (|3.77p 
can be formulated as an upper bound for vi. Note that the pressure as a 
unified characteristics of all Euclidean Gibbs states makes senses for translation 
invariant models only. Thus, the notions mentioned in Definition 13.71 are not 
applicable to the versions of the model which do not possess this property. 
The main result of this subsection is contained in the following statement. 

Theorem 3.22 Let the model be as just described. Let also the condition Ij3.71\ l 
with defined by the equation \3. 77p and J defined by \3. 76'[ ) be satisfied. Then 
for every (3 > /?*, where /3* is defined by the equation \3. 71\ ), the model has 
a phase transition in the sense of Definition \3.1\ If the model is translation 
invariant, the long range order and the first order phase transition take place at 
such p. 

Proof: The proof is made by comparing the model under consideration with a 
reference model, which is the scalar model with the nearest neighbor interaction 
of intensity (|3.76p and with the anharmonic potential (|3.78p . Thanks to the 
condition (|3.77p . the reference model is more stable; hence, the phase transition 
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in this model implies the same for the model considered. The comparison is 
conducted by means of the correlation inequalities. 

The reference model is translation invariant and hence can be defined by its 
local periodic Hamiltonians 

where for a box A, E is the same as in (|3.32j) : iJ^ ar is as in For this 

model, we have the infrared estimate (|3.55|) with v = 1. Let us obtain the lower 
bound, see (|3.58[) . To this end we use the inequalities p.59[) . (|3.60|) and obtain 

r 

< a + 2&«+y2s(2s-l)fc( s VM0)] 2(s ~ 1) ) low (3.82) 

s=2 

s=2 V '■ 

Here v l ™ is the periodic Gibbs measure for the model (|3.81[) . To get the second 
line we used the Gaussian upper bound inequality, see page 1031 in [54] and 
page 1372 in [Hj, which is possible since all b^ s \ s > 2 are nonnegative. The 
solution of the latter inequality is 

(MO)) 2 ) , ow >tf*. (3.83) 

Then the proof of the phase transitions in the model (|3.81|) goes along the line 
of arguments used in proving Theorem 13.201 Thus, for (3 > /?*, (wefi))^™ > 0, 

where /x+ w is the corresponding maximal Euclidean Gibbs measure, see Propo- 
sition [HEl But, 

M0)) M+ > M0)) M i r , (3.84) 

see Lemma 7.7 in [53]. At the same time (oty(0)) M = for any periodic /i € G t ', 
which yields the result to be proven. □ 

3.5 Phase transition in the scalar model with asymmetric 
potential 

The phase transitions proven so far have a common feature - the spontaneous 
symmetry breaking. This means that the symmetry, e.g., rotation invariancc, 
possessed by the model and hence by the unique element of Q t is no longer 
possessed by the multiple Gibbs measures appearing as its result. In this sub- 
section, we show that the translation invariant scalar version o the model (|l.ljl , 
(|1.2p has a phase transition without symmetry breaking. However, we restrict 
ourselves to the case of first order phase transitions, see Definition 13.71 The 
reason for this can be explained as follows. The fact that D^ e , does not decay 
to zero as \l — l'\ — > +oo, see p.!8[) . implies that ^ is non-ergodic only if is 
symmetric. Otherwise, to show that /i is non-ergodic one should prove that the 
difference D^, — (fe)^ • (fe 1 )^ does not decay to zero, which cannot be done by 
means of our methods based on the infrared estimate. 

In what follows, we consider the translation invariant scalar version of the 
model (jl.ip . (|1.2[) with the nearest neighbor interaction. The only condition 



(MO)) 2 ) 
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imposed on the anharmonic potential is (|2.2[) . Obviously, we have to include 
the external field, that is the anharmonic potential is now V(u) — hu. Since 
we are not going to impose any conditions on the odd part of V, we cannot 
apply the GKS inequalities, see [3 Ell, the comparison methods are based on, 
see p.84|) . In view of this fact we suppose that the interaction is of nearest 
neighbor type. Thus, for a box A, the periodic local Hamiltonian of the model 
has the form (|3~5Tj) . 

In accordance with Definition 13. 71 our goal is to show that the model pa- 
rameters (except for h) and the inverse temperature (3 can be chosen in such a 
way that the set TZ, defined by (|2 . 100(1 . is non-void. The main idea on how to 
do this can be explained as follows. First we find a condition, independent of 
h, under which Df,, does not decay to zero for a certain periodic fi. Next we 
prove the following 

Lemma 3.23 There exist h±, h- < h + , which may depend on the model pa- 
rameters and (3, such that the magnetization iS.lOl]) has the property: 

M(h)<0, for freft c n(-oo,/i_); M(h) > 0, for h G TZ C n (h+ + oo). 

Thereby, if TZ were void, one would find h* G (h-,h+) such that M(/i*) = 0. 
At such h*, the aforementioned property of would yield the non-ergodicity 
of [i and hence the first order phase transition, see Theorem 13.221 

In view of Corollary [231 D^, does not decay to zero if (|3.16|) holds with big 
enough By Proposition 13.181 the lower bound (|3.16|l can be obtained from 
the estimate (|3.61|) . The only problem with the latter estimate is that it holds 
for h = 0. 

Lemma 3.24 For every (3 > and 9, there exist positive m* and J», which 
may depend on (3 > and 9 hut are independent of h, such that, for any box A 
and any /i£i, 

( MO)] 2 ) /cr >9, if J > J* and m > m*. (3.85) 
Proof: For ft 6 1, we set 



\ h (duj) = —f-expih u(T)dT A(dw), (3.86) 




Np = I exp ( h I lo(t)&t ) X(duj) 

where A is as in (|2.102[) . Then for ±h > 0, we get the estimate p.64j) in the 
following form 

0Jd{ [lu { {0)} 2 > i/P e, > Jde 2 + log X h [^(^c)] , (3.87) 

where B ± (e, c) is as in (|2.108|> . (|2.109|) . Let us show now that, for ±h > 0, 

X h [B ± (e, c)] > A [B± (e, c)] . (3.88) 

For h > 0, let I(ui) be the indicator function of the set Ct(n;c), see (|2.106p . 
For S > and t G R, we set 

t < c, 

= ^ ( f - c)/ 5 * G (c,c + (5], 
c > c + <5. 
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Thereby, 



h{to) 



dcf 



\{ L 5 Hk/3/n)} 



k=0 



By Lebesgue's dominated convergence theorem 



Np\ h 



C+(n;c) 



C/3 



lim 



no j Cl 



I{lj) exp [ h J cj(r)drj X(doj) 

Is(uj)exp h / w(r)dr A(dw). 



(3.89) 



As the function Ig is continuous and increasing, by the FKG inequality, see 
Theorem 6.1 in [3], it follows that 

f I s (x)ex-p ( h / u(r)dT ] A(dw) > JV^ / I ff (w)A(dw). 

jc^ \ io y Jc'/n 

Passing here to the limit we obtain from (|3.89[) 

C+(n;c)] > \[c+(n;c) 

which obviously yields (|3.88p . For h < 0, one just changes the signs of h and u>. 
Thereby, we can rewrite p.87p as follows, c.f., (|3.64p . 

(N(0)] 2 >„ r >e 2 + [\o gl {m)]/(5Jd. 

Then one applies the arguments from the very end of the proof of Lemma 13.161 
□ 

Proof of Lemma \3.23t Suppose that h > 0. Then restricting the integration in 
Until to [B+(e,c)] A , we get 



Pa 



p A c \h) > hps + log N P + ^ £ 2 Y<Jm>+ log \{B + (e,c)} (3.90) 



'6A 



> h{3e + logiV/3 + log7(m). 



As the right-hand side of the latter estimate is independent of A, it can be 
extended to the limiting pressure p{h). For any positive h 6 TZ C , by the convexity 
of p(h) one has 

M(h) > [p(h)-p(0)]/(3h 

> £+^{-p(0) + logAf / 3 + log 7 (m)}. 

Picking big enough h we get the positivity stated. The negativity can be proven 
in the same way. □ 

Now we are at a position to prove the main statement of this subsection. 

Theorem 3.25 Let the model be scalar, translation invariant, and with the 
nearest-neighbor interaction. Let also d > 3. Then for every (3, there exist 
to* > and J* > such that, for all m > to* and J > J*, there exists h* G R, 
possibly dependent on m, (3, and J, such that p'(h) gets discontinuous at h*, 
i.e., the model has a first order phase transition. 
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Proof: Let m* be as in (|2.110[) and J* , 6 be as in Lemma^OH Fix any (3 > and 
m > m*. Then, for J > J*, the estimate p. 851) holds, which yields the validity 
of (|3.74[) for all boxes A with such /?, m, and v — 1. Thereby, we increase J, if 
necessary, up to the value at which (|3.75p holds. Afterwards, all the parameters, 
except for h, are set fixed. In this case, there exists a periodic state fj, € G t such 
that the first summand in (|3.15p is positive; hence, D^, does not decay to zero 
as \l — i'\ — > +oo, see (|3.12l) and (|3.15p . If p(h) is everywhere differentiable, 
i.e., if 1Z = 0, then by Lemma T3. 231 there exists h* such that M(h*) = 0; hence, 
the state \x with such h* is non-ergodic, which yields \Q t \ > 1 and hence a first 
order phase transition. Otherwise, 1Z ^ 0. □ 

3.6 Comments 

• Subsection \3.1l According to Definition 13. 1[ the phase transition corre- 
sponds to the existence of multiple equilibrium phases at the same values 
of the model parameters and temperature. This is a standard definition 
for theories, which employ Gibbs states, see [30] . In the translation invari- 
ant case, a way of proving phase transitions can be to show the existence 
of non-ergodic elements of C*'. For classical lattice systems, it was re- 
alized in [55] by means of infrared estimates. More or less at the same 
time, an alternative rigorous theory of phase transitions in classical lat- 
tice spin models based on contour estimates has been proposed. This is 
the Pirogov-Sinai theory elaborated in [HT], see also [70]. Later on, this 
theory was essentially extended and generalized into an abstract sophisti- 
cated method, applicable also to classical (but not quantum) models with 
unbounded spins, see [M] and the references therein. 

For quantum lattice models, the theory of phase transitions has essential 
peculiarities, which distinguish it from the corresponding theory of classi- 
cal systems. Most of the results in this domain were obtained by means of 
quantum versions of the method of infrared estimates. The first publica- 
tion in which such estimates were applied to quantum spin models seems 
to be the article [Mj ■ After certain modifications this method was applied 
to a number of models with unbounded Hamiltonians [5] [15 ] IT6 ] [231144], 160]. 
In our approach, the quantum crystal is described as a system of 'classi- 
cal' infinite dimensional spins. This allows for applying here the original 
version of the method of infrared estimates elaborated in [35] adapted 
to the infinite dimensional case, which has been realized in the present 
work. Among others, the adaptation consists in employing such tools as 
the Garsia-Rodemich-Rumsey lemma, see [29] . This our approach is more 
effective and transparent than the one used in [5J [THl EH HI] . It also allows 
for comparing the conditions (|3.16[) . (|3.17p with the stability conditions 
obtained in the next section. 

In the physical literature, there exist definitions of phase transitions alter- 
native to Definition 13.11 based directly on the thermodynamic properties 
of the system. These are the definition employing the differentiability 
of the pressure (Definition 13. 7\ which is applicable to translation invari- 
ant models only), and the definition based on the long range order. The 
relationship between the latter two notions is established by means of 
the Griffiths theorem, Proposition 13.91 the proof of which can be found 
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in [24] . For translation invariant models with bounded interaction, non- 
diffcrcntiability of the pressure corresponds to the non-uniqueness of the 
Gibbs states, see [33 We failed to prove this for our model. 

In the language of limit theorems of probability theory, the appearance of 
the long range order corresponds to the fact that a new law of large num- 
bers comes to power, see Theorem 13.91 and the discussion preceding Defi- 
nition [nHU The critical point of the model corresponds to the case where 
the law of large numbers still holds in its original form (in the translation 
invariant case this means absence of the first order phase transitions), but 
the central limit theorem holds true with an abnormal normalization. For 
a hierarchical version of the model (|l.ip . (|1.2p . the critical point was de- 
scribed in [3] . Algebras of abnormal fluctuation operators were studied in 
[20j . In application to quantum crystals, such operators were discussed 
in [80l [81] , where the reader can find a more detailed discussion of this 
subject as well as the corresponding bibliography. 

• Subsection VS.'Bi As was mentioned above, the method of infrared esti- 
mates was originated in |28j . The version employed here is close to the 
one presented in [5B]. We note that, in accordance with the conditions 
(|3.9p . (|3.16|) . and ()3.17|) . the infrared bound was obtained for the Duhamel 
function, see (|3.55p . rather than for 

e'eA 

which was used in J5[ [151 HE1 S3] • 

• Subsection \3.3\ The lower bound (|3.58p was obtained in the spirit of [531 
R)U] . The estimate stated in Lemma r3.16l is completely new; the key element 
of its proving is the estimate (|2.105p . obtained by means of Proposition 
12.31 The sufficient condition for the phase transition obtained in Theorem 
13.201 is also new. Its significant feature is the appearance of a universal 
parameter responsible for the phase transition, which includes the particle 
mass m, the anharmonicity parameter and the interaction strength J. 
This is the parameter on the left-hand side of (|3.7ip . The same very 
parameter will describe the stability of the model studied in the next 
section. Theorem 13.211 is also new. 

• Subsection \3.4\ Here we mostly repeat the corresponding results of [54] . 
announced in [53] . 

• Subsection \3.5[ The main characteristic feature of the scalar model stud- 
ied in [5[ CD3 QH [23 El [SO], as well the the one described by Theorem 
13.221 was the ^-symmetry broken by the phase transition. This sym- 
metry allowed for obtaining estimates like (|3.83p . crucial for the method. 
However, in classical models, for proving phase transitions by means of the 
infrared estimates, symmetry was not especially important, see Theorem 
3.5 in [55] and the discussion preceding this theorem. There might be two 
explanations of such a discrepancy: (a) the symmetry was the key element 
but only of the methods employed therein, and, like in the classical case, 
its lack does not imply the lack of phase transitions; (b) the symmetry 
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is crucial in view of e.g. quantum effects, which stabilize the system, see 
the next section. So far, there has been no possibility to check which 
of these explanations is true. Theorem 13.251 solves this dilemma in favor 
of explanation (a). Its main element is again an estimate, obtained by 
means of the Garsia-Rodemich-Rumsey lemma. The corresponding result 
was announced in [39] . 

4 Quantum Stabilization 

In physical substances containing light quantum particles moving in multi- welled 
potential fields phase transitions are experimentally suppressed by application 
of strong hydrostatic pressure, which makes the wells closer to each other and 
increases the tunneling of the particles. The same effect is achieved by replacing 
the particles with the ones having smaller mass. The aim of this section is to 
obtain a description of such effects in the framework of the theory developed 
here and to compare it with the theory of phase transitions presented in the 
previous section. 

4.1 The stability of quantum crystals 

Let us look at the scalar harmonic version of the model Ijl.lj) - a quantum 
harmonic crystal. For this model, the one-particle Hamiltonian includes first 
two terms of (|1.2j) only. Its spectrum consists of the eigenvalues = (n + 
l/2)y / a/m, neffo. The parameter a > is the oscillator rigidity. For reasons, 
which become clear in a while, we consider the following gap parameter 

A^ = mm(£jr - E^). (4.1) 

nGN 

Then 

A h ™ = a = mAl aT . (4.2) 

The set of tempered Euclidean Gibbs measures of the harmonic crystal can be 
constructed similarly as it was done in section [21 but with one exception. Such 
measures exist only under the stability condition (|2.71[) . which might now be 
rewritten 

Jo < mAl m . (4.3) 

In this case, Q t is a singleton at all f3, that readily follows from Theorem 12.141 
As the right-hand side of (|4.3[) is independent of to, this stability condition is 
applicable also to the classical harmonic crystal which is obtained in the classical 
limit to — > +oo, see [§]. According to (|2.2|) the anharmonic potentials Vt have a 
super-quadratic growth due to which the tempered Euclidean Gibbs measures of 
anharmonic crystals exist for all Jo . In this case, the instability of the crystal is 
connected with phase transitions. A sufficient condition for some of the models 
described in the previous section to have a phase transition may be derived from 
the equation (|3.75p . It is 

2/?J^/(/3/4toi9») > J(d), (4.4) 
which in the classical limit to — > +oo takes the form 



> J{d). 
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The latter condition can be satisfied by picking big enough j3. Therefore, the 
classical anharmonic crystals always have phase transitions - no matter how 
small is the interaction intensity. For finite m, the left-hand side of (|4.4[) is 
bounded by 8m?9* J, and the bound is achieved in the limit (3 — > +oo. If for- 
given values of the interaction parameter J, the mass m, and the parameter 
which characterizes the anharmonic potential, this bound does not exceed 
J(d), the condition (|4.4|) will never be satisfied. Although this condition is only 
sufficient, one might expect that the phase transition can be eliminated at all 
j3 if the compound parameter &md 2 J is small enough. Such an effect, if really 
exists, could be called quantum stabilization since it is principally impossible in 
the classical analog of the model. 



4.2 Quantum rigidity 

In the harmonic case, big values of the rigidity a ensure the stability. In this 
subsection, we introduce and stugy quantum rigidity, which plays a similar role 
in the anharmonic case 

Above the sufficient condition ([4.4)1 for a phase transition to occur was ob- 
tained for a simplified version of the model ([1.1)1 , (jl.2p - nearest neighbor inter- 
actions, polynomial anharmonic potentials of special kind ([3.78p . ect. Then the 
results were extended to more general models via correlation inequalities. Like- 
wise here, we start with a simple scalar version of the one-particle Hamiltonian 
HEm , which we take in the form 

H m = -^f + \? + V(q), (4.5) 
where the anharmonic potential is, c.f., ()3.78[) . 

V(q) = b^q 2 + b^q 4 + ■ ■ ■ + b^q 2r , 6< r > > 0, rGN\{l}. (4.6) 



The subscript m in (|4.5| indicates the dependence of the Hamiltonian on the 
mass. Recall that H m acts in the physical Hilbert space L 2 (K). Its relevant 
properties are summarized in the following 

Proposition 4.1 The Hamiltonian H m is essentially self-adjoint on the set 
(R) of infinitely differentiate functions with compact support. The spectrum 
of H m has the following properties: (a) it consists of eigenvalues E n , n G No 
only; (b) to each E n there corresponds exactly one eigenfunction tp n £ L 2 (R); 
(c) there exists 7 > 1 such that 

n" 1 E n — > +00, as n — > +00. (4-7) 

Proof: The essential self-adjointness of H m follows from the Sears theorem, see 
Theorem 1.1, page 50 of [T7] or Theorem X.29 of [64) . The spectral properties 
follow from Theorem 3.1, page 57 (claim (a)) and Proposition 3.3, page 65 (claim 
(b)), both taken from the book [17]. To prove claim (c) we employ a classical 
formula, see equation (7.7.4), page 151 of the book [77], which in our context 
reads 

-V2^ f " y/E n - V(u) du = n + i + O ( - ) , (4.8) 
tt Jo 2 \nj 
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where n, and hence E n , are big enough so that the equation 

V(u) = E n 

have the unique positive solution u n . Then 

7T 



,r+l 







H) 









(4.9) 



(4.10) 



where 



ITtV - + a /2)* 2 ^V" 1 )^" 1 ). 



Note that n (l) = 1 for all n, which follows from 

E„ 



Thus, 



1, 



as n-t +oo. 



Thereby, we have 



del 



^(j> n {t)-t 2r dt 

Then combining (|4~T^|) with ([43]) and (|4~TT|) we get 

"^r(| + i) 



V 7 ! - t 2r dt 



(4.11) 



(4.12) 



-Em 



r 6 « - 


l/(r+l) 


_{2mY_ 





r(l)r(i) 



0(1) 



(4.13) 



which readily yields (|4.7|) with any 76 (1, 2r/(r + 1)). □ 

Thus, in view of the property (|4.13[) we introduce the gap parameter 

A m = min(-E n - 25 n _i), 



(4.14) 



and thereby, c.f., ([472]) . 

ft m = mZ\^, (4.15) 

which can be called quantum rigidity of the oscillator. One might expect that 
the stability condition for quantum anharmonic crystals, at least for their scalar 
versions with the anharmonic potentials independent of £, is similar to (|4.3|) . 
That is, it has the form 

Jo<n m . (4.16) 



4.3 Properties of quantum rigidity 

Below / ~ g means that lim(//g) = 1. 

Theorem 4.2 For every r£ N, the gap parameter A m , and hence the quan- 
tum rigidity lZ m corresponding to the Hamiltonian \4--5\j , J^.6| ), are continuous 
functions of m. Furthermore, 

A m ~ A Q m- r ^ r+l \ Tl m ~ A 2 m-^-^ r+1 \ m -> 0, (4.17) 

with a certain Aq > 0. 
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Proof: Given a > 0, let U a : L 2 (M.) — > L 2 (R) be the following unitary operator 

(J7 Q V0 (x) = y/ai()(ax). (4-18) 

Then by HQ) 



Fix any mo > and set p — (m/roo) 1 '' r " 1 " 1 ', a = p 1 / 2 . Then 



H m = U- x H m U a = p- r T(p), (4.19) 



where 



T(p) = H mo +Q(p) (4.20) 
1 -p 2 + p r -\b^ + a/2)q 2 + p r - 2 b^q 4 + ■■• + b {r) q 2r , 



2mo 

Q(p) = (p-i)[ Pr ^(p)(b^ +a/2)q 2 (4.21) 
+ Pr - 2 (p)b^q 4 + ■ ■ ■+p r -s(p)b^q 2s + ■■■ + b^-Vq 2 ^' , 



and 



p k (p) = 1 + p + p z + ■■■ + p k -\ (4.22 



As the operators H m , H m , are unitary equivalent, their gap parameters (|4.14p 
coincide. The operators H m and T(p), p > possess the properties established 
by Proposition ^. II In particular, they have the property (|4.7|) with one and the 
same 7. Therefore, there exist e > and k £ N such that for \p — 1| < e, the 
gap parameters (|4.14|) for if m and T(p) are defined by the first k eigenvalues 
of these operators. As an essentially self-adjoint operator, T(p) possesses a 
unique self-adjoint extension T(p), the eigenvalues of which coincide with those 
of T(p). Furthermore, for complex p, T(p) is a closed operator, its domain 
Dom[T(p)] does not depend on p. For every %j> £ Dom[T(p)], the map C 3 
C h-> f(C,)ip £ L 2 {R) is holomorphic. Therefore, {f (p) | |p - 1| < s} is a self- 
adjoint holomorphic family. Hence, the eigenvalues 0„(,o), n £ No of T(p) are 
continuous functions of p £ (1 — e, 1 + e), see Chapter VII, §3 in the book [4*0] , 
At p = 1 they coincide with those of H mo . Since we have given k £ N such that, 
for all p £ (1 — e, 1+ e), 

mm [6„(p) - e„_i(»] = min [9„(p) - Q n -l(p)] , 
neN n£{l,2,...,fe} 

the function 

Z(p) d ^ f minp-'- [0 n (p) - e n _i(p)] (4.23) 
is continuous. But by (|4. 19[) 

Z\ m = A ((m/m ) 1/(r+1) ) , (4.24) 

which proves the continuity stated since mo > has been chosen arbitrarily. 
To prove the second part of the theorem we rewrite (|4.20p as follows 

T(P)=H£1+R(p), (4.25) 
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where 



1 



= —v 2 + b^a 2r 
m ° 2m P + q ' 

and 

R(p) = p + a /2)q 2 + p r -%Wq A + ■■■ + fcC'" V (r_1) ) . 

Repeating the above perturbation arguments one concludes that the self-adjoint 
family {T(p) \ \p\ < e} is holomorphic at zero; hence, the gap parameter of (|4.25j) 
tends, as p — > 0, to that of Hml, i.e., to Aq. Thereby, the asymptotics (|4.17p 
for A m follows from (|4.19p and the unitary equivalence of H m and H m . □ 

Our second result in this domain is the quasi-classical analysis of the param- 
eters (|4.14j) . (|4.15p . Here we shall suppose that the anharmonic potential V has 
the form (j4~6|) with > for all s = 2, . . . , r - 1, c.f., (|3J8| . We remind that 
in this case the parameter > is the unique solution of the equation (|3.79|) . 

Theorem 4.3 Let V be as in {3. 78\j . Then the gap parameter A m and the 
quantum rigidity lZ m of the Hamiltonian \4-5\ ) with such V obey the estimates 

A - * 2^? U - ± 4^f (4 ' 26) 



Proof: Let g m be the local Gibbs state (|2.5|) corresponding to the Hamiltonian 
(|4.5|) . Then by means of the inequality (|3.59|) and the Gaussian upper bound 
we get, see (I3.82|) . 

a + 2b& + t> {g m (q 2 )) >0, 

by which 

Qm(q 2 )>#*. (4.27) 

Let tp n , n G No be the eigenfunctions of the Hamiltonian H m corresponding to 
the eigenvalues E n . By Proposition 14. 1[ to each E n there corresponds exactly 
one ip n - Set 

Qnn> = (^n,#n')L2(R), n,Tl' G N . 

Obviously, Q n n — for any n G No. Consider 

T(r, r') = g m [ 9 exp (-(r' - r)H m ) gexp (-(r - T r )H m )\ , r, r' G [0, (3], 

which is the Matsubara function corresponding to the state g m and the operators 
Fi=F 2 = q. Set 

u(k)= r(0,r)cosfcTdr, fc G /C = {(2tt//3)k |k G Z}. (4.28) 
Jo 

Then 

n,n'=0 

x {exp(-/3# n ') - exp(-f3E n )} , 
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where Z m = trace cxp(—/3H m ). The term (E n — E n /) 2 in the denominator can 
be estimated by means of (|4.14[) , which yields 

1 1 + °° 

^ ^ u-TT'T- E IQwl 2 {E n -E n .) (4.30) 
x {exp(-(3E n ) - exp(-f3E n ,)} 

i 



By this estimate we get 

g m (q 2 ) = r(0,0) = |2«(fc) (4.31) 

Combining the latter estimate with (|4.27[) we arrive at 
Z\ m tanh(/3Z\ m /2) < l/(2m0.), 

which yields l|4.26p in the limit (3 — > +oo. □ 

Now let us analyze the quantum stability condition (|4.16[) in the light of 
the latter results. The first conclusion is that, unlike to the case of harmonic 
oscillators, this condition can be satisfied for all Jo by letting the mass be small 
enough. For the nearest-neighbor interaction, one has Jo = 2dJ; hence, if (|4.16[) 
holds, then 

8dm&lJ < 1. (4.32) 
This can be compared with the estimate 

Zdm-dlJ > dj{d), (4.33) 



guaranteeing a phase transition, which one derives from 

For finite d, dj{d) > 1, see Proposition 13. 19) hence, there is a gap between 
the latter estimate and (|4.32[) . which however diminishes as d — > +oo since 

lira dj(d) = 1. 

d — > + oo 

In the remaining part of this section, we show that for the quantum crystals, 
both scalar and vector, a stability condition like (|4.16p yields a sufficient decay 
of the pair correlation function. In the scalar case, this decay guaranties the 
uniqueness of tempered Euclidean Gibbs measures. However, in the vector case 
it yields a weaker result - suppression of the long range order and of the phase 
transitions of any order in the sense of Definition 13.71 The discrepancy arises 
from the fact that the uniqueness criteria based on the FKG inequalities are 
applicable to scalar models only. 
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4.4 Decay of correlations in the scalar case 

In this subsection, we consider the model (|1.2[) which is (a) translation 

invariant; (b) scalar; (c) the anharmonic potential is V(q) = v(q 2 ) with v being 
convex on R + . 

Let A be the box (|2.63p and A* be its conjugate l|3.2p . For this A, let 

K^(r,r') ^ (MrU'ir'))^ (4.34) 

be the periodic correlation function. Recall that the periodic interaction poten- 
tial Jffif was defined by (|2.78p . For the one-particle Hamiltonian (|1.2|) . let u(k) 
be as in P~2"8"]) . 

Theorem 4.4 Let the model be as just describes. If 

m(0)J o < 1, (4.35) 

then 

K A r i\ ^ 1 e^[t(p,e-£')+tk(T-T')} 

w/iere 

J A = E 4', r A (p) = J A - ]T 4, exp[i(p, * - 1% (4.37) 
reA reA 

Proof: Along with the periodic local Gibbs measure (|2.82j) we introduce 

^ Cr (cMi) (4.38) 

= at^Ta exp \ \ E J w{ui,u t ,) h i - [ J2v(uj e (T))dT 1 XA(dw A ), 

JV A [t) ( f,i'£A ^ ^£A J 

where t £ [0, 1] and N^° r (t) is the corresponding normalization factor. Thereby, 
we set 

X u ,(Ty\t)=(ut(T)w i ,(/)) vrm , £,£'£A. (4.39) 
By direct calculation 

^-X u ,(T,r'\t) (4.40) 
at 



l r 

o E J ^2 / ^«'^ 2 (7-,r',T",T"|f)dr" 
E J ^ f X Ul {r,T"\t)X l2l ,{r"y\t)dr". / 



where 



Re 1 i 2 e 3 e 4 ,(n ) T 2 ,T 3) T4\t) = (^(ti)^^)^^)^^)) ^« ( .| t) 

- (o;^(Ti)^ 2 (r2))„p« ( .| t) • (^ 3 (r3)w <4 (r4)) 1 ^«(.|t) 

- (^ 1 (n)^ 3 (r 3 )) 1/ p»(.| t) ■ (^ 2 (T 2 )^ 4 (r 4 ))^ ( .| t) 

- (^ 1 (n)^ 4 (r4)) 1/ P«(.|t) ■ (^ 2 (T2)^ 3 (T3))^«(.| t ). 
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By the Lebowitz inequality, see [6], we have 



^2/3/4^1, T 2 ,T3,r 4 |t) < 0, (4.41) 



holding for all values of the arguments. Let us consider (|4.40j) as an integro- 
differential equation subject to the initial condition 



X U '(t,t'\0) = 5 u >T(t,t') = (<5«7/3)^u(fc)cosfc(T-T'). (4.42) 
Besides, we also have 

X u ,(t,t'\1) = K a ,(t,t'\ P ). (4.43) 

Along with the Cauchy problem (|4.40p . (|4.42H let us consider the following 
equation 

d 



Of 



Y U >(T,T>\t)= £ 



26A 



T 1 
J I 



|A| 



Y Ul (T,r"\t)Y e2i ,(r",r'\t)dT", (4.44) 



where e > is a parameter, subject to the initial condition 
Y^(t,t'|0)=J£«/(t,t'|Q) 

= (foV/3)$^«(fc)cOSfc(T-T / ). 



(4.45) 



Let us show that under the condition ()4.35|) there exists eo > such that, for 
all e S [0,£o), the problem (|4.44p . (|4.45[1 . t e [0, 1], has the unique solution 



Y, 



ex P [l(p,£-f) +lk(T-T')} 



P\M tt f r - t{J£ + s6 p , ] + tr A (p) ■ 



(4.46) 



where Jo, T A (p) are the same as in (|4.37|) and 5 Pj o is the Kronecker symbol with 
respect to each of the components of p. By means of the Fourier transformation 

Y u ,(r,r'\t) = -L J2 X) ? (f' fc W ex PWP^-0+^(r-r')](4.47) 



Y(p, k\t) 



V / r«,(r,r'|i)exp[- l (p,£-f)-zfc(r-r')]dr', 



we bring l|4.44p . (|4.45p into the following form 

^Y(p, k\t) = [j A (p) + eS pfi ] ■ [?(p, fc|t)] 2 , ?(p, fc|0) = u(k), (4.48) 



where, see (|4.37|) . 

j a (p) = 4' cx P [»(p, * - 1')} = J A - r A ( P ). 



(4.49) 



'6A 



Clearly, J A < J , \J A {p)\ < J A , and u{k) < fi(0). Then in view of (|435]) . one 
finds £o > such that, for all e € (0, £0)1 the following holds 



J A {p) + e5 Pi0 u{k) < 1, 
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for all p G A* and k E K.. Thus, the problem (|4.48|) can be solved explicitly, 
which via the transformation (|4.47[) yields (|4.46p . 
Given 9 G (0, 1), we set 



Y e { e V(T,T'\t)=Y u ,(T,T'\t + 6), te[O,l-0]. (4.50) 

Obviously, the latter function obeys the equation (|4.44|) on t e [0, 1 — 8] with 
the initial condition 

Y$(t,t'\0) = Y u ,(t,t'\0) > Y u ,{t,t'\Q) = X u ,(t,t'\0). (4.51) 

The latter inequality is due to the positivity of both sides of (|4.44|) . Therefore, 

Y^(r,r'\t)>0, (4.52) 

for all £, £' G A, t, t' G [0, (3], and t G [0, 1 - 9]. 

Let us show now that under the condition (|4.35[) , for all 8 G (0,1) and 

£ G (0,£ ), 

X w {T,r'\t) < Y$\ T ,T'\t), (4.53) 
also for all £, £' G A, r, r' G [0, 0\, and i G [0, 1 - 9). To this end we introduce 

Zf t ,{T,r'\t) d ^ Y^( T ,T'\t)±X u ,(T,T'\t), te [0,1-0]. (4.54) 
Then one has from ([4~40ll . (|4T4]) 



: Z u ,{T,T'\t) (4.55) 



5 E J tij {Z+jT,T"\t)Z^(T\T"\t) 



£i,£ 2 eA 

+Z tti (r,r"|t)Z+ 2 (r',r"|t)}dr" 
rP 



where Sw(t, r'\t) stands for the first term on the right-hand side of (|4.40p . By 
(1431) and 63U) 

Z«, (t,t'\0) = F«,(r,r'|0) - X«/(r,r'|0) > 0, (4.56) 

which holds for all £,£' G A, t,t' G [0,/3]. For every £,£' G A, both F«'(r,r'|i), 
Xh'(t, r'\t) and, hence, Z%,(t, r'|t) are continuous functions of their arguments. 
Set 

C(i) =inf{Z«,(r,r'|i) K,/G A, r, r' G [0, /?]} . (4.57) 

By (|4.56p . it follows that £(0) > 0. Suppose now that C(^o) = at some 
t G [0, 1 - 9} and > for all t G [0, t ). Then by the continuity of Z^i, 
there exist £,£' G A and r, r' G [0, /3] such that 

Z~,(t, r%) = and Z^,(t, r'\t) > for all t < t . 

For these £,£' G A and r, r' G [0,/3], the derivative (d/dt)Z~,(T, r'|f) at t = i 
is positive since on the right-hand side of (|4.55[) the third term is positive and 
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the remaining terms are non-negative. But a differentiable function, which is 
positive at t € [0,to) and zero at t = to, cannot increase at t = to- Thus, 
C(t) > for all t £ [0, 1 - 6], which yields (|4.53p . By the latter estimate, we 
have 

X u ,(t,t'\1-9)<Y u ,(t,t'\1) 

1 exp[i(p,£-£')+ik(T-T')] 

All the function above depend on 9 and e continuously. Hence, passing here to 
the limit 9 = e \ and taking into account (14.43)) we obtain (|4.36[) . □ 

By means of Proposition 12.211 the result just proven can be extended to all 
periodic elements of G t - For /j, £ Q t , we set 

K^{t,t') = {u 1 {t)^ i {t')) ii . (4.58) 



Theorem 4.5 Lei i/ie stability condition |^.i6| ) 6e satisfied. Then for every 
periodic fi 6 £/ , f/ie correlation function J^.5ff| j /ias i/ie bound 

K^(t,t') < Y u >(tS) (4.59) 
def 1 /" expMp.l-lQ+tfcCT-T')] 

Wj^c •/(-*.*]" [«(fc)]- 1 -io+ r(p) p ' 

w/iere 

r(p)=J -^J«/expKp^-O], pe(-^,7r] d . (4.60) 

£' 

T/ie same bound has also the correlation function K^,(t,t'), where fiQ G (? is 
f/ie same as m Provosition \2. 1§1 

Remark 4.6 By gTggp , [w^)]" 1 > m(zi^ + fc 2 ). T/ie -upper 6ownd in UJty 
with [■u(fc)]~ 1 ] replaced by m([Z\ har ] 2 + fc 2 ) turns into the infinite volume corre- 
lation function for the quantum harmonic crystal discussed at the beginning of 
subsection \4-l\ Thus, under the condition ^4-35] ) the decay of the correlation 
functions in the periodic states is not less than it is in the stable quantum har- 
monic crystal. As we shall see in the next subsection, such a decay stabilizes 
also anharmonic ones. 

For T(p) ~ 2"o |p| 2 , To > 0, as p — > 0, the asymptotics of the bound in (|4.59p 
as y/\£ — £'\ 2 + \t — t'| 2 — * +oo will be the same as for the d + 1-dimensional 
free field, which is well known, see claim (c) of Proposition 7.2.1, page 162 of 
[3"T] . Thus, we have the following 

Proposition 4.7 If the function j4.60\ ) is such that T(p) ~ To|p| 2 , To > 0, as 
p — > 0, the upper bound in has an exponential spacial decay. 
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4.5 Decay of correlations in the vector case 

In the vector case, the eigenvalues of the Hamiltonian (|4. 5|) are no longer simple; 
hence, the parameter (|4.14p definitely equals zero. Therefore, one has to pick 
another parameter, which can describe the quantum rigidity in this case. If the 
model is rotation invariant, its dimensionality v is just a parameter. Thus, one 
can compare the stability of such a model with the stability of the model with 
v = 1. This approach was developed in see also [HI [35]. Here we present 
the most general result in this domain, which is then used to study the quantum 
stabilization in the vector case. 

We begin by introducing the corresponding class of functions. A function 
/ : R — > R is called polynomially bounded if f(x)/(l + \x\ k ) is bounded for some 
k £ N. Let J- be the set of continuous polynomially bounded / : R — > R which 
are either odd and increasing or even and positive. 

Proposition 4.8 Suppose that the model is rotation invariant and for alii G A, 
AgL, Vi(x) = vg(\x\ 2 ) with vt being convex on R + . Then for any t\, . . . , r n € 

[o,/3], 4, •••,4 e A, j = /i,.../„ eF, 

(/i(a4f(ri))---/„(4f (T n )))„ A < (hMri))---fnMr n ))h A , (4.61) 

where D\ is the Euclidean Gibbs measure &2.33\) of the scalar model with the same 
Jui as the model considered and with the anharmonic potentials Vi(q) — vi(q 2 ). 

By this statement one immediately gets the following fact. 

Theorem 4.9 Let the model be translation invariant and such as in Proposition 
\4-S\ Let also A m be the gap parameter J^. 11$ of the scalar model with the same 
interaction intensities Ju> and with the anharmonic potentials V{q) = v{q 2 ). 
Then if the stability condition {^.16\j is satisfied, the longitudinal correlation 
function 

K^{ry) = {Jf?{r)Jf^)) ll , j — 1,2, ... (4.62) 

corresponding to any of the periodic states n € Q , as well as to any of the 
accumulation points of the family {7TA(-|0)}AgL, obeys the estimate |^.5ff[ ) in 
which u{k) is calculated according to fr4-29\ ) for the one- dimensional anharmonic 
oscillator of mass m and the anharmonic potential v(q 2 ). 

4.6 Suppression of phase transitions 

From the 'physical' point of view, the decay of correlations (|4.59[) already cor- 
responds to the lack of any phase transition. However, in the mathematical 
theory, one should show this as a mathematical fact basing on the definition of 
a phase transition. The most general one is Definition 13.11 according to which 
the suppression of phase transitions corresponds to the uniqueness of tempered 
Euclidean Gibbs states. Properties like the differentiability of the pressure, c.f., 
Definition 13. 7\ or the lack of the order parameter, see Definition 13. 5[ may also 
indicate the suppression of phase transitions, but in a weaker sense. The aim of 
this section is to demonstrate that the decay of correlations caused by the quan- 
tum stabilization yields the two-times differentiability of the pressure, which in 
the scalar case yields the uniqueness. This result is then extended to the models 
which are not necessarily translation invariant. 
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In the scalar case, the most general result is the following statement, see 
Theorem 3.13 in [SI] , 

Theorem 4.10 Let the anharmonic potentials Ve be even and such that there 
exists a convex function v : R+ — > R, such that, for any Ve, 

Ve(xg) — v(x t ) < Ve(xe) — v(x 2 ) whenever x\ < x\. (4.63) 

For such v, let A m be the gap parameter of the one-particle Hamiltonian 1 1.1]) 
with the anharmonic potential v(q 2 ). Then the set of tempered Euclidean Gibbs 
measures of this model is a singleton if the stability condition J4.-?6'| ] involving 
A m and the interaction parameter Jq of this model is satisfied. 

The proof of this theorem is conducted by comparing the model with the trans- 
lation invariant reference model with the anharmonicity potential V(q) = v(q 2 ). 
By Proposition 12.18) for the model considered and the reference model, there 
exist maximal elements, /x+ and /i™ f , respectively. By means of the symmetry 
Ve(q) = Vi(— q) and the FKG inequality, one proves that, for both models, the 
uniqueness occurs if 

{wt(0))tf = 0, M0)) M+ = 0, for all I. (4.64) 

By the GKS inequalities, the condition (|4.63p implies 

0<M0)W < M0))^ f , (4.65) 

which means that the reference model is less stable with respect to the phase 
transitions than the initial model. The reference model is translation invariant. 
By means of a technique employing this fact, one proves that the decay of 
correlations in the reference model which occurs under the stability condition 
(|4.16j) yields, see Theorem Ol 

Mo))^- = o, 

and therefrom ()4.64[) by (|4.65|) . The details can be found in [54] , 

As was mentioned above, in the vector case we did not manage to prove that 
the decay of correlations implies the uniqueness. The main reason for this is 
that the proof of Theorem 14.101 was based on the FKG inequality, which can 
be proven for scalar models only. In the vector case, we get a weaker result, by 
which the decay of correlations yields the normality of thermal fluctuations. To 
this end we introduce the fluctuation operators 

Q^-L^gCa A <g L, j = l,...,v. (4.66) 
V l A l eeA 

Such operators correspond to normal fluctuations. 

Definition 4.11 The fluctuations of the displacements of oscillators are called 
normal if the Matsubara functions h2.12]) for the operators F\ = Q^ 1 ' , ■ ■ ■ , F n = 
Q(i™) ; remain bounded as A f L. 
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If A is a box, the parameter (|3.28[) can be written 

^ = iyrr rft B) (3) (T,r')drdr'. (4.67) 

Thus, if the fluctuations are normal, phase transitions of the second order (and 
all the more of the first order) do not occur. 

Like in the proof of Theorem 14. 9[ the model is compared with the scalar 
ferromagnetic model with the same mass and the anharmonic potential v(q 2 ). 
Then the gap parameter A m is the one calculated for the latter model. 

Theorem 4.12 Let the model be the same as in Theorem \4-^\ and let the sta- 
bility condition involving the interaction parameter Jq of the model and the gap 
parameter A m corresponding to its scalar analog be satisfied. Then the fluctua- 
tions of the displacements of the oscillators remain normal at all temperatures. 

4.7 Comments 

• Subsection \4-l\ In an ionic crystal, the ions usually form massive com- 
plexes the dynamics of which determine the physical properties of the 
crystal, including its instability with respect to structural phase transi- 
tions, see |21j . Such massive complexes can be considered as classical 
particles; hence, the phase transitions are described in the framework of 
classical statistical mechanics. At the same time, in a number of ionic 
crystals containing localized light ions certain aspects of the phase tran- 
sitions are apparently unusual from the point of view of classical physics. 
Their presence can only be explained in a quantum-mechanical context, 
which points out on the essential role of the light ions. This influence of 
the quantum effects on the phase transition was detected experimentally 
already in the early 1970's. Here we mention the data presented in [TSUTo] 
on the KDP-type ferroelectrics and in [S7J on the YBaCuO-type super- 
conductors. These data were then used for justifying the corresponding 
theoretical models and tools of their study. On a theoretical level, the 
influence of quantum effects on the structural phase transitions in ionic 
crystals was first discussed in the paper [65j , where the particle mass was 
chosen as the only parameter responsible for these effects. The conclu- 
sion, obtained there was that the long range order, see Definition 13. 51 gets 
impossible at all temperatures if the mass is sufficiently small. Later on, 
a number of rigorous studies of quantum effects inspired by this result as 
well as by the corresponding experimental data have appeared, see [SHI 152"] 
and the references therein. Like in [65] . in these works the reduced mass 
(|1.3[) was the only parameter responsible for the effects. The result ob- 
tained was that the long range order is suppressed at all temperatures in 
the light mass limit m — > 0. Based on the study of the quantum crystals 
performed in [21 E] [5j [7J [9] , a mechanism of quantum effects leading to the 
stabilization against phase transitions was proposed, see [5]. 

• Subsection \4-ty According to [S] the key parameter responsible for the 
quantum stabilization is lZ m — mA^, see (|4.15[) . In the harmonic case, 
mAf n is merely the oscillator rigidity and the stability of the crystal cor- 
responds to large values of this quantity. That is why the parameter 
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toZ\^ was called quantum rigidity and the effect was called quantum sta- 
bilization. If the tunneling between the wells gets more intensive (closer 
minima), or if the mass diminishes, mA^ n gets bigger and the particle 'for- 
gets' about the details of the potential energy in the vicinity of the origin 
(including instability) and oscillates as if its equilibrium at zero is stable, 
like in the harmonic case. 

• Subsection \4-3[ Theorems l4.2l and l4.3l are new. Preliminary results of this 
kind were obtained in [3J 150] . 

• Subsection \4-4\ Theorems 14.41 Wl>[ 14.71 were proven in [45j . 

• Subsection ^. 5\ Various scalar domination estimates were obtained in [471 

• Subsection \4- 6] ' Theorem 14.101 was proven in [54] . The proof of Theorem 
14.121 was done in . The suppression of abnormal fluctuations in the 
hierarchical version of the model (jl.ip , (|1.2p was proven in [2] . 
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